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ABSTRACT 


) 


The  development  of  ship  hull  forms  of  decreased  wave-making 
resistance  by  mathematical  methods  has  received  some  investigation  in 
the  past  few  years  but  little  application*  However,  it  appears  that 
much  more  application  can  be  obtained  by  use  of  improved  mathematical 
methods  recently  developed.  The  theory  and  effectiveness  of  work  done 
by  other  authors  is  reviewed.  In  addition,  the  mathematical  machinery 
necessary  to  calculate  hull  forms  by  a  method  of  steep  descent  is 
developed  and  applied  to  a  simple  example.  The  necessary  resistance 
equations  are  formulated.  A  closed  form  solution  is  obtained  for  the 
smoothly  varying  portion  of  the  wave-making  resistance  of  an  assemblage 
of  sources  and  sinks  traveling  near  the  surface  of  the  water.  An 
Integral  for  the  interference  terms  which  comprise  the  fluctuating 
portion  of  the  wave -making  resistance  is  also  presented  in  the  shape  of 
a  Laplace  transform.  It  is  concluded  that  much  more  extensive  use  may 
be  made  of  mathematical  methods  to  improve  hull  forms  than  has  been  the 
case  heretofore. 


iv 


I.  Introduction 


The  purpose  of  this  paper  Is  to  examine  what  can  be  done  to  calculate  ^ 
the  shapes  of  ship  hulls  which  will  have  less  wave-making  resistance  than 
those  developed  by  ordinary  drafting  techniques.  The  work  done  In  this 
field  will  be  reviewed  and  It  will  be  shown  to  contain  evidence  that  the 
mathematical  approach  has  a  large  potential  for  making  improvements.  This 
has  not,  however,  been  exploited  for  two  reasons.  First,  the  mathematical 
techniques  used  so  far  are  not  flexible  enough  to  be  aoplled  to  many 
practical  cases.  Second,  the  wave-making  resistance  formula  on  which  these 
techniques  are  based  has  not  produced  precise  predictions  of  experimental 
results,  and  this  has  cast  doubt  on  the  validity  of  computations  derived 

from  it.  It  will  be  pointed  out  in  this  paper  that  a  way  of  bringing  theory 

and  experiment  into  agreement  which  has  been  worked  out  by  Inui  for  certain 
special  cases  can  probably  be  extended  for  more  general  application.  In 
addition,  a  mathematical  method  of  more  general  applicability  than  those 
used  until  now  will  be  discussed  --  the  method  of  steep  descent  --  and  the 
mathematical  machinery  and  reformulation  of  the  wave -making  resistance  equations 
necessary  to  permit  its  use  will  be  outlined.  At  this  point  it  will  be  conclu¬ 
ded  that  the  way  is  now  open  to  more  extensive  use  of  mathematical  techniques 

for  the  calculation  of  practical  ship  forms  of  decreased  wave-making  resistance. 

A  few  words  on  the  history  of  the  wave -making  resistance  problem  may 
serve  to  explain  these  comments,  and  show  how  the  stage  has  been  set  for  an 
advance.  In  the  past,  although  there  has  been  a  great  deal  of  experimental 
measurement  of  the  wave-making  resistance  of  ships,  there  has  been  little  use 
of  mathematical  methods  to  find  ways  to  decrease  this  resistance.  Since 
William  Froude  demonstrated  by  his  experiments  the  general  correctness  of  the 
modeling  law  for  wave-making  resistance  In  the  1870's,  it  has  been  possible 
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for  naval  architects  to  predict  from  model  tests  the  wave-making  resistance  of 
the  ships  corresponding  to  the  particular  models.  This  prediction  added  to  an 
estimate  of  the  "frictional  resistance"  has  been  accurate  enough  to  establish  ' 
the  horsepower  required  to  drive  the  full-sized  ships.  Many  thousands  of 
models  have  been  towed  and  good  predictions  made.  Naval  architects  have  also 
understood  the  relationship  between  the  overall  dimensions  of  a  ship  and  its 
wave-making  resistance.  It  has  been  recognized,  for  example,  that  if  a  given 
amount  of  displacement  is  placed  in  a  long  ship,  this  ship  will  probably  have 
less  wave-making  resistance  than  a  shorter  ship  of  the  same  displacement.  Until 
recently,  however,  there  has  been  little  application  of  detailed  understanding 
of  the  theory  of  wave  production  by  ships'  hulls  to  the  improvement  of  tnose 
hulls.  So  little  systematic  improvement  was  made  over  many  years  that  a 
certain  set  of  hull  designs  developed  early  in  the  twentieth  century,  Taylor's 
Standard  Series,  was  long  taken  as  a  standard  of  goodness.  If  a  hull  design 
had  as  little  wave  resistance  as  Taylor's  models,  it  was  considered  a  very 
good  design.  Most  designs  had  more. 

Despite  the  lack  of  progress  in  application  of  theory  to  Improving  hull 
designs,  a  considerable  amount  of  theoretical  understanding  of  wave  resistance 
was  developed  by  a  few  investigators,  starting  with  J.  Mlchell,  who  published 
a  classic  paper  on  the  subject  in  1898  [l].  Mlchell  derived  an  expression  for 
the  wave -making  resistance  of  a  thin  hull  moving  on  the  surface  of  an  ideal 
fluid.  Sir  Thomas  Havelock  followed  with  a  long  series  of  papers  in  which 
this  basic  theory  was  applied,  expanded,  and  reformulated  in  a  more  tractable 
form,  which  has  been  shown  [2]  to  be  equivalent  to  the  original  theory  of 
Mlchell.  He  considered  the  hull  to  be  the  set  of  closed  streamlines  generated 
by  a  set  of  moving  sources  and  sinks.  Unfortunately,  two  factors  caused  naval 
architects  to  disregard  the  Mlchell -Have lock  theory.  First,  the  results  of 
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experiment,  although  qualitatively  the  same  as  the  predictions  of  the  theory, 
never  agreed  with  it  exactly.  Second,  the  application  of  the  theory  to 
practical  hull  forms  required  a  great  deal  of  calculation.  > 

Several  developments  in  the  past  ten  years  appear  likely  to  permit  far 
more  use  of  theoretical  analysis  In  the  development  of  hull  forms  of  less  wave 
resistance.  Not  the  least  important  is  the  general  availability  of  high  speed 
computers.  In  addition  to  this,  however,  Takao  Inui  in  Japan  has  laid  the 
groundwork  for  a  more  accurate  calculation  of  the  actual  wave  resistance. 

First,  he  has  found  that  by  describing  the  hull  by  the  streamlines  generated 
by  the  set  of  sources  and  sinks  employed  in  Havelock's  theory  rather  than 
using  a  simpler  approximation  to  the  hull,  he  can  get  much  better  agreement 
with  experiment  than  before.  Second,  by  using  a  small  number  of  semi -empirical 
parameters  to  account  for  the  sheltering  effect  of  the  hull,  the  motion  of  the 
wake,  and  the  effect  of  linearizing  approximations,  he  can  get  precise  agree¬ 
ment  between  experiment  and  theory  [3}  .  In  addition,  he  has  demonstrated 
theoretically  and  experimentally  that  it  is  possible  in  the  case  of  certain 
hull  shapes  to  add  a  spherical  appendage  at  bow  and  stern  in  such  a  location 
as  nearly  to  cancel  the  bow  waves  and  stern  waves  produced  by  the  hull  M. 

In  this  country  Weinblura  developed  tables  which  have  permitted  the  use  of 
Rltz's  method  to  find  improvements  in  hulls  which  can  be  described  by  a  limited 
set  of  polynomials  C^3  *  Even  though  the  calculated  resistance  does  not  agree 
exactly  with  the  experimental  resistance,  several  hulls  developed  in  this  way 
have  been  found  experimentally  to  be  considerably  better  than  those  of  the 
same  dimensions  developed  in  the  ordinary  way. 

Using  a  different  approach,  Karp,  Kotik,  and  Lurye  have  succeeded  in 
applying  the  classical  calculus  of  variations  to  the  problem  of  finding  the 
strut  of  minimum  wave-making  resistance  •  The  particular  strut  is  one  which 
can  be  described  by  a  distribution  of  dipoles  along  a  plane  of  finite  length 
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but  Infinite  depth.  the  distribution  being  Independent  of  the  depth  end 
described  between  the  ends  of  the  plane  by  functions  lntegrable  when  multi¬ 
plied  by  e  certain  Green's  function.  While  the  results  are  of  limited  practical 
importance.  they  do  demonstrate  that  a  solution  does  exist  in  this  case  —  a 
matter  concerning  which  there  had  been  some  question. 

Despite  these  encouraging  developments  there  has  not  yet  been  developed 
a  general  approach  which  will  permit  the  systematic  improvement  of  practical 
hull  forms.  It  is  not  always  possible  .  for  example  .  to  adopt  a  form  which 
can  be  expressed  as  a  polynomial  .  or  to  Install  a  sphere  beneath  the  bow. 

In  any  case  most  of  the  shape  of  the  ship  will  be  dictated  by  what  it  has  to 
carry,  stability,  and  other  considerations,  so  that  features  to  reduce  wave 
resistance  must  be  made  compatible.  There  seems  to  be  no  reason  why  an  approach 
cannot  be  developed  which  will  do  this.  This  paper  will  first  review  the  work 
which  has  been  done  toward  finding  ways  to  reduce  the  wave-making  resistance 
of  hulls,  and  will  then  suggest  how  a  more  general  approach,  using  a  method 
of  steep  descent,  may  be  applied  to  improve  any  hull  form  whose  major 
characteristics  have  already  been  set. 


4 


II.  Derivation  and  Discussion  of  Wave  Resistance  Equations 
Both  theory  and  experiment  show  wave-making  resistance  of  a  ship  results 
from  the  system  of  waves  which  the  ship  leaves  behind  as  it  travels  over  the  ^ 
surface  of  the  water.  Although  at  a  constant  ship  speed  the  waves  made  by 
the  ship  form  a  pattern  which  appears  to  move  with  the  ship,  this  pattern 
actually  extends  farther  and  farther  aft  with  time. 

The  energy  required  to  increase  its  size  is  provided  by  the  ship.  Since 
the  equations  which  describe  both  the  pattern  and  the  energy  required  to  pro^ 
duce  it  have  been  developed  by  others,  and  the  development  is  quite  difficult, 
only  as  much  of  the  derivation  will  be  included  here  as  to  make  clear  the 
limitations  in  their  use. 

1.  ASSUMPTIONS  AND  EQUATIONS  OF  MOTION 
The  assumptions  in  the  first  part  of  the  following  discussion  are  those  £  7 J 
employed  by  Lunde.  They  will  be  modified  later  in  certain  cases  to  account 
for  some  semi -empirical  corrections  made  by  Inui.  The  notation  is  similar  to 
that  used  by  Lunde,  although  not  exactly  the  same.  It  will  be  used  through¬ 
out  this  paper  for  the  sake  of  consistency,  even  though  it  differs  consider¬ 
ably  from  that  of  many  of  the  authors  referred  to.  The  assumptions  are  these: 

a.  The  fluid  is  incompressible,  homogeneous,  and  of  zero  viscosity. 

b.  The  motion  is  irrotational .  It  can  therefore  be  described  by  a  velocity 
potential  0  such  that  the  fluid  velocity  vector  is  given  by 

q  ■  •  V(l.  Hence  ifq*iu+jv  +  fcw,  where  i,  j,  and  k  are  unit  vectors 
parallel  to  the  cartesian  coordinates  x,  y,  z,  the  fluid  velocity  has  the 
components 

u  -  -  0X,  v  -  -  0y,  w  »  -  0Z. 

c.  The  wave  height  is  small  in  comparison  with  the  wave  length,  so  that  the 
wave  slope  is  a  small  quantity. 


d.  The  ship  has  been  moving  in  a  straight  line  on  calm,  infinitely  deep  water 
for  an  infinitely  long  time. 

a.  The  motion  set  up  by  the  ship  can  be  approximated  by  the  motion  produced  ^ 
by  a  set  of  sources  and  sinks  or  doublets. 

f.  The  wave-making  resistance  is  Independent  of  the  frictional  and  eddy 
resistance  of  the  ship  and  can  be  calculated  separately. 

The  notation  will  employ  cartesian  coordinates  throughout.  The  following 
conventions  will  be  employed: 

The  axis  of  x  will  be  in  the  mean  surface  of  the  water,  and  it  will  be 
oriented  so  that  x  is  positive  in  the  direction  in  which  the  ship  is  moving, 
y  is  positive  to  starboard,  and  z  is  positive  upward.  The  origin  of  the 
x-coordinates  is  taken  in  the  ship  and  the  x-axis  moves  with  the  ship. 

The  steady  speed  of  the  ship  in  the  x  direction  will  be  c. 

The  acceleration  of  gravity  will  be  g. 

The  parameter  KQ  -  g/c*. 

The  elevation  of  the  surface  of  the  water  will  be  and  this  will  be 
positive  in  the  up  direction.  It  will  take  the  value  zero  at  z  -  0  (the  mean 
position  of  the  surface). 

For  differential  equations  we  have  the  following: 

aH  +g2g  +a2g  -  o 

Syl  jrg  ...(2.1) 

This  is  the  equation  of  continuity. 

3^0  +  KJ3J|  "  0  on  z*  0. 

3?  %  ...(2.2) 

This  provides  the  equation  of  wave  motion  on  the  surface. 

In  addition,  it  is  possible  to  relate  the  elevation  of  the  surface  to  the 
velocity  potential  by  the  following  equation: 


2.  THE  VELOCITY  POTENTIAL  OF  A  MOVING  SOURCE 

Using  the  assumptions  discussed  above ,  an  expression  for  the  velocity 
potential  has  been  derived  both  by  Peters  and  Stoker  and  by  Lunde,  using  ^ 
different  approaches*  In  this  section  their  results  will  be  shown  to  be 
equivalent* 

It  is  necessary  that  the  expression  for  0  be  a  solution  of  both  (2.1) 
and  (2.2),  that  it  approach  cx  at  infinity,  that  it  give  a  pattern  which  in 
the  vicinity  of  the  ship  moves  with  the  ship,  and  that  the  water  be  undisturbed 
at  a  relatively  small  distance  ahead  of  the  ship.  All  these  requirements  are 
met  by  solving  (2.1)  and  (2.2)  as  an  initial  value  problem  and  then  letting 
the  time  go  to  infinity.  Peters  and  Stoker  used  this  device  of  letting  the  time 
go  to  infinity  rather  than  the  technique  of  adding  an  artificial  frictional 
resistance  used  by  Lunde.  In  this  way  Peters  and  Stoker  obtained  a  Green's 
function  which  can  be  interpreted  as  the  velocity  potential  of  an  isolated 
moving  source  [,8j.  Let  m  be  the  strength  of  the  source,  and  let  (h ,  k,  -f)  be 
its  location  in  cartesian  coordinates.  Then  if  we  convert  the  Green's  function 
into  Lunde's  notation  we  get  the  following  expression  for  the  velocity  potential 


0  ■ 

m  -  _ _ m  4* 

^/(x-h)^  +  (y-k)2  +  (z+f)2  y(x-h)2  +“(y.k)2  +  (z-f)2 

+  4mg  Re  S  eK(z-f)  eiK(x-h)cos9  cos  [k(y-k)sinOj  dK  dQ  « 

^  °  L  g  -  K  cos^G 


m  -  m 

y(x-li)2,  T  (y-k)2  +  (z+i :)2  jJx-W)7  +  (y-k)2  +  (z-f)2 

-  furtKo  Re  ?  ^  eK(z-f)  e-iK(x-h)cos9  cos  [K(y-k)3infl]  aec20  dK  dO 
*  o  l  K  -  Kosec20 


...(2.4) 
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The  contour  L  goes  from  the  origin  to  infinity  along  the  real  axis  but  is 
deformed  above  the  real  axis  in  the  vicinity  of  K  *  Kosec20. 

In  the  case  of  a  ship  there  will  always  be  symmetry  about  the  x-z  plane/ 
so  that  for  each  source  m  at  (h,  k,  -f)  there  will  be  another  at  (h,  -k,  -f ) . 

It  follows  that  the  velocity  potential  will  also  be  symmetrical  in  the  x-z 
plane  so  that 

(x,0,z)  «  0  ...(2.5) 

ay 

From  this  it  follows  that  if  we  describe  a  ship  by  a  set  of  sources  distri¬ 
buted  symmetrically  with  respect  to  the  x-z  plane  we  will  get  a  velocity 
potential  corresponding  to  any  symmetrical  pair  which  is  the  sum  of  two 
velocity  potentials.  This  is  simply 

0*  -  0(x,y,z;  h,k,-f)  +  0(x,y,z;  h,-k,-f)  ...(2.6) 

If  we  perform  the  addition  indicated  in  equation  (2.6)  we  obtain  the 
velocity  potential  explicitly. 

0*  « 

_ m _  -  _ m  __  + 

>/  (x-h)^  +  (y-k)2  +  (z+f)2  J  (x-h)^  +  (y-k)^  +  (z-f)2 

+  _ m _  - _ m _ _ 

J  (x-h)2  +  "(y+k)2  +  (z+f)2  >/  (x-h)2  +  (y+k)2  +(z-f)2  ~~ 

n/2 

“  Re  \  J  eK(z-f)  o-iR(x-h)cos0^2)co3(Kysin9)coa (KksinQ)  sec2fl  dK  dO 
**  o  L  K-KqScc^O 

...(2.7) 

The  contour  integration  of  (2.7)  gives  for  the  last  term  the  value  2^*  Res(Kosec^0) 
Hence  if  we  call  the  denominators  of  the  first  four  terms  of  (2.7)  ,  r2 ,  r ^ , 

and  r^,  we  have  the  following  results: 


8 


9*  -  -  (  1  -  l  +1  -1  )  - 
ri  r2  r3  r4 

-  16 ^  cKo(*-f)s«c  ®s£n^Ko(x»|l)secQj  cosjVoysinQsec2©}  x  N 

o 

x  cos  £lC0k  sinQsec2©}  sec29  dO  .*.(2.8) 

This  is  the  wave  pattern  produced  by  a  symmetrical  pair  of  sources ,  one  to 
port  and  one  to  starboard  of  the  centerline  of  a  ship.  We  observe  that  at 
s  *  0  the  terms  in  r^  through  r^  cancel  out  and  leave  us  with  the  integral 
term  only. 

We  may  compare  this  with  Lunde's  expression  for  the  elevation  of  the 
surface  at  a  great  distance  aft  of  the  ship.  To  do  this  we  take 

X-  -(c/g )9$*  -  16  m  K„2(c/g)  \  eKf,^Z’f^8ec2®co8rK0(x-h)sece  '|  x 

'  5 X  0  L 

x  cos  £  K0ysin9sec2Qj  cos£lC0k  sinO  sec2f)J  sec^O  dO 

...(2.9) 


For  k  ■  0  this  Is  equivalent  to  his  expression  (7.10): 

I2  C  mgCosT  KQ(x-h8)  sec  ol  cos  T  K0(y-k8)sec2©  sinoT  x 
'  c  -*372  s-1 

x  exp(-K  f  sec^O)  sec^O  dO 

...(2.9a) 

The  depth  -z  will  be  zero  in  both  equations.  The  change  in  range  of  Integra* 
tion  from  -fl/2  to  7F/2  down  to  0  to  T/2  together  with  the  fact  that  equation 
(2.9)  assumes  two  sources  —  one  port  and  one  starboard  --  each  of  strength  m, 
while  Lunde's  equation  refers  to  only  one  source,  explains  the  difference  of  a 
factor  of  4  between  the  equations.  The  identity 
K02(c/g)  -  K0(g/c2)(c/g)  «  Ko/c 

completes  the  demonstration. 
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Sine*  Peters  end  Stoker  obtained  their  result  without  recourse  to  artificial 
friction  levs ,  the  fact  that  their  equation  is  the  sene  as  Lunds' s  shows 
that  Lunde's  later  conclusions  following  this  equation  do  not  depend  on  the  ^ 
existence  of  such  friction  in  water*. 

3.  THE  RESISTANCE  INTEGRAL 

If  a  ship  travels  into  quiet  water  it  will  cause  an  Increase  In  the  wave 
energy  in  the  region  into  which  it  advances.  By  calculating  the  rate  of  change 
of  the  wave  energy  in  a  large  region  containing  the  ship  together  with  the 
eaount  of  wave  energy  which  crosses  the  boundary  of  the  region  it  is  possible 
to  calculate  the  anount  of  energy  which  the  ship  converts  into  waves.  Since 
we  know  the  notion  of  the  water  everywhere  once  we  have  the  velocity  potential, 
it  is  certainly  possible  to  calculate  the  energy  put  into  the  water  by  the 
ship  in  this  wanner. 

Let  there  be  a  fixed  large  area  of  the  water  surface  into  which  the 
ship  is  advancing.  It  is  sufficient  to  bound  this  by  two  infinite  planes , 
one  well  forward  of  the  ship  and  one  well  aft  of  it,  each  at  right  angles  to 
the  direction  of  the  ship's  motion.  We  will  call  these  respectively  plane 
A  and  plane  B.  If  we  let  E(A)  and  E(B)  be  the  rate  at  which  wave  energy  crosses 
into  the  area  across  boundaries  A  and  B,  W(A)  and  W(B)  be  the  rate  at  which 
work  is  done  on  the  fluid  within  these  boundaries  at  the  boundaries ,  R  the 
wave  resistance  and  c  the  speed  of  the  ship,  then  we  can  write  the  energy 
balance  on  the  large  region: 

*ln  discussing  such  comparisons  as  these,  it  is  worth  remenbering  that  the 
definitions  of  m  used  by  some  authors  (Inui,  for  example)  differ  from  the  one 
used  by  Lunde  by  a  factor  4fl/c.  This  results  in  a  corresponding  factor  c/4? T 
in  the  coefficient  of  the  expression  for  the  wave  elevation  . 
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R  c  -  K(A)  -  B(B)  +  W(A)  +  W(B)  -  0  ...(2.10) 

Bine*  the  plane  far  ahead  of  the  ahip  la  undisturbed  by  th*  ahlp'a  notion, 
w  can  aat  E(A)  “  W(A)  ■  0.  We  than  can  write  th*  wav*  raaiatanca  of  the  , 
ahip  as 

R-  (1/c)  [b(B)  -W(B)J  ...(2.11) 

By  substituting  for  B(B)  and  W(B)  their  equivalent  in  terns  of  th*  velocity 
potential  and  going  through  some  extensive  nanipulatlon,  Lunde  obtains  finally 
an  sxprssslon  for  the  wave  resistance -^7^. 


R  -  167T*K02  )  (Pe2  +  P02  +  Q.2  +  Q02)  aec3«  dO 

*  o 

-Kof**c2®  , 

Pe«  $  6e  cos (KQh  secO)  cos  (KQk  sinO  aec*9)  ds 

s 

-Kgfsec2® 

Pe-  ^6*  sin(Rohaec0)sin(Koksin0  sec2®)  ds 

s 

-Kofsec20 

$6e  sln(Rohaec0)cos(Koksln0  sec2®)  ds 

a 

-ILfaec2® 

Qcm  ^  6  e  cos(Kohaec0)sln(Kokain0  sec2®)  ds 

a 


...(2.12) 


...(2.13) 


Here  6  is  the  source  density  at  any  point  (h,  k,  -f)  on  the  hull  or  within 
its  boundary  and  s  indicates  integration  over  the  volume  of  the  hull.  The 
quantity  p  is  the  density  of  water. 

It  la  possible  to  simplify  this  expression  even  further  in  the  case  that 
the  hull  is  symmetrical  about  the  centerline  plane.  Then  the  following  result 
can  be  used: 


16yi 


?<H 


■f  Q2)  tec*  9  <19 


cot 

tin 


J^Ko(hcoa0  +  ksln®)sec2®J  exp(-KQfsec20)  ds 


...(2.14) 

...(2.13) 


11 


It  should  bo  emphasized  that  in  this  case  it  is  necessary  to  carry  out  the 
integration  over  both  halves  of  the  ship,  since  the  simplification  has  made 
the  result  correct  only  when  the  two  halves  are  added.  This  should  be 
apparent  from  an  examination  of  equation  (2.15),  which  is  not  ea  even 
function  of  k. 

Another  relationship  which  makes  the  interactions  more  obvious  can  be 
derived  from  equations  (2.12)  and  (2.13).  Here  we  will  substitute  a  finite 
sum  for  the  integral. 

6  (*)  - »  £  «r 

r-1 

If  we  substitute  this  in  equations  (2.12)  and  (2.13),  combine  terms  and 
simplify,  we  finally  get  the  following  expression: 


16*J»K02  )  £  m  2  ^  e"2K°fr8ec2®  sec3®  d®  + 

(  r-1  o 


+  2  £  £  mrmB  S  e“Ko(fr+fs)*ec20  coaf K0(hc-hs)seco'] 

s-r+1  r-1  o  L  J 

jf K^kj-kg)  sinOsec2®  J  sec3®  d®  ^ 


x  cos 


...(2.16) 


Equation  (2.16)  may  be  put  in  the  form  of  a  Laplace  transform.  He  will 
let  t+1  ■  sec2®.  Then  equation  (2.16  can  be  written  in  the  following  form: 


R  -  16i^>K02 


'IVr*  (l+t)1/2  t“1/2dt 


li‘'Yrrb 

*  I  Pr*  j  .  Pt*  (l+t)1/2  t*1/2  CO.Jqc,<lM)1,^CO.Jij.tl/2a+t)1/2|llA 


...(2.17) 


Here  we  have  set 

®rr  "  2Kofr5  Prs  *  K0(£r+ft):  qr,  ■  Mhr‘ha>;  q£g  “  K0<kr“ks^ ' 
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It  should  bo  oboorvod  that  tho  expression  is  on  ovon  function  of  k9  and  there- 
foro  it  provides  the,  same  vsluo  for  the  Interaction  between  a  source  on  the 
centerline  end  either  of  two  symmetric  sources  off-center  on  either  side* 
Equations  (2.14)  and  (2*15)  provide  a  different  value  for  the  Interaction  * 
tern  between  a  source  on  the  centerline  and  a  source  on  the  port  side  than 
they  do  for  the  Interaction  term  between  the  source  on  the  centerline  and 
the  symmetrical  source  on  the  starboard  side,  and  provide  a  correct  result 
only  when  all  the  interaction  terms  are  summed.  This  complicates  interpreta¬ 
tion.  On  the  other  hand,  equations  (2.16)  and  (2.17),  which  are  actually  more 
general,  provide  terms  which  can  be  Interpreted  directly  as  interactions 
between  the  sources  m^  and  m#.  In  consequence,  equation  (2.17)  will  be  used 
in  much  of  the  discussion  later. 

The  first  Integral  of  equation  (2.17)  can  be  evaluated  in  closed  form. 

Let  us  call  the  resistance  corresponding  to  this  first  term  R^1)  ,  and  that 
corresponding  to  the  second  term  and  let  R  be  their  sum: 

R  -  R(1)  +  R<2>  ...(2.18) 

We  aay  now  write  r(D  in  closed  font.  We  will  let  rO)  be  described  by  a  sun. 

R(1)  -  £  R  ...(2.19) 

r«l 

Rrr  -  prr  1«r2®’3Pr,f/2kl(Prr/2)  ...(2.20) 

The  function  k^(prr/2)  is  Bateman's  function.  Derivation  of  equation  (2.20) 
is  outlined  in  appendix  1.  For  very  small  KQ,  which  corresponds  to  very  large 
c,  the  value  of  Rrr  is  small.  As  KQ  increases,  Rrr  Increases,  reaches  a 
maximum,  and  then  decreases.  The  term  rM,  which  is  the  smoothly  varying 
part  of  the  wave-making  resistance,  is  a  function  only  of  the  depth  of  the 
sources,  not  of  their  position  along  the  hull  or  their  distance  outboard  of  the 
centerline • 
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There  seems  to  be  no  closed  fora  velue  for  the  teres  which  comprise 
R<2>.  On  the  other  hand,  they  ere  closely  related  In  fora  to  the  terns  which 
coeprlse  R^  .  Let  us  write  R^2^  as  s  sue.  > 

n  n»l 

R<2>  -  2  £  £Rre  .  ...(2.21) 

s*r+l  r»l 


Then  we  can  write 


•r.  •  “r'oS-.f 

x  CO.  TqJ,  t1/J(l«)lrt  J  <t  ...(2.22) 

It  is  possible  to  show  If  q^#  ■  0  that  if  the  integral  is  taken  over  a  finite 
interval  from  the  origin  to  any  zero  of  cos  k.  the  error  will 

be  leF^  than  the  value  of  the  Integral  from  that  zero  to  the  next  one.  The 
same  can  be  shown  if  -  0,  q'  j  0,  except  that  the  zeros  referred  to  are 
those  of  cos  ^qrit1^2(l+t)1^2J  .  The  demonstration  of  this  is  shown  in 
appendix  1*  When  both  qra  and  q^fl  differ  from  zero9  the  estimation  of  the 
error  Incurred  by  terminating  the  integration  is  less  simple  but  It  is  possible 
to  get  an  upper  limit  by  observing  that  for  all  t?t0>  l<(l+t)*/2t-l/2  <(i+l/to)*^* 
Consequently  the  error  will  be  less  than 


J  s-Prat(l+l/t0)1/2  dt|  -  J (l/pr#)  e"Pr,t°(l+l/to)1/2 

Also  it  is  apparent  from  examination  of  equations  (2,17)  through  (2,22)  that 
the  magnitude  of  the  interaction  term  between  any  pair  of  sources  mr  and  mg 
must  be  less  than  or  equal  to  the  magnitude  of  the  sum  of  the  resistance  of 
mr  and  m§  taken  separately* 


2R 


rs 
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<  R  +  R 

*  rr  as 


...(2.23) 


The  equality  sign  can  hold  only  if  the  two  aources  are  superposed.  It  is9 
of  course ,  possible  to  have  a  negative  value  for  Rr##  The  inequality  shows 
that  no  pair  of  sources  can  Interact  with  each  other  in  such  fashion  as  to  ^ 
produce  a  negative  total  resistance. 

4.  EMPIRICAL  CORRECTIONS 

Comparison  of  the  wave-making  resistance  as  calculated  by  equation  (2.14) 
with  the  resistance  as  measured  in  model  basins  has  produced  good  qualitative 
agreement,  but  in  general  has  not  produced  exact  quantitative  agreement*  In 
particular  it  has  been  found  that  the  humps  and  hollows  of  the  curve  of  resist¬ 
ance  vs.  speed,  although  of  approximately  the  right  separation  in  speed  are 
somewhat  displaced  from  their  proper  position,  and  the  calculated  humps  and 
hollows  appear  exaggerated  when  compared  with  the  experimental  ones.  In  addi¬ 
tion,  there  have  been  some  systematic  differences  in  the  magnitude  of  the 
resistance.  In  view  of  the  approximations  which  have  been  made  in  deriving 
the  resistance  equations  (no  viscosity,  squares  of  velocities  other  than 
ship  speed  negligible,  small  slope  of  waves,  and  so  forth)  these  differences 
are  not  altogether  surprising.  Fortunately,  an  explanation  of  the  major 
differences  has  been  produced  by  Takao  Inui.  In  the  case  of  a  number  of 
models  whose  source  distributions  could  be  described  by  simple  continuous 
functions  he  has  found  it  possible  to  correct  for  these  differences  with  a 
relatively  small  number  of  semi -empirical  parameters  [3j.  It  appears  possible 
to  extend  these  techniques  to  other  distributions. 

Inui1 s  first  discovery  was  that  the  agreement  between  theory  and  experi¬ 
ment  could  be  significantly  improved  if  the  hull  form  were  found  by  plotting 
the  closed  streamlines  produced  by  the  set  of  sources  and  sinks  employed  to 
describe  the  hull,  rather  than  using  the  approximation  employed  previously.  The 
approximation  used  by  earlier  investigators  was  this: 
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6  -  -  c_  dk 

2r  *h  •••(2.24) 

Here  $  is  the  source  density ,  end  h  end  k  ere  the  x  end  y  coordinates  of  the 

> 

hull  boundary*  This  approximation  is  only  satisfactory  for  very  thin  ships , 
but  Inul  showed  j^3j  that  it  is  unsatisfactory  where  the  beam  is  as  great  as 
one -twentieth  the  length  --  a  ship  still  thinner  thaq  any  used  in  practice* 

An  example  of  extreme  differences  in  results  between  using  approximation  (2.24) 
and  the  actual  plotting  of  the  streamlines  can  be  obtained  by  an  illustration 
of  the  results  for  a  two-dimensional  source  forward  of  a  two-dimensional  sink* 

In  the  approximation  (2*24)  this  amounts  to  an  infinitely  deep  rectangular 
barge;  however,  if  the  streamlines  are  traced  a  somewhat  longer,  more  fair 
shape  is  produced  --  still,  of  course,  in  the  form  of  an  infinitely  deep  strut* 


Source 

motion 

Source  (*Hn) 


Two-dimensional  source  and  sink  by  Mlchell's  approximation  (2.24)* 


His  next  discovery  [3  ]  was  that  in  the  case  of  hulls  which  have  both  a 
well-defined  bow  wave  and  a  well-defined  stern  wave,  with  no  waves  originating 
between  bow  and  stern,  he  could  bring  the  calculated  and  observed  humps  in 
the  resistance  curve  together*  For  purposes  of  calculation  of  the  fluctuating  part 
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of  the  resistance  he  essuaed  that  the  stern  wave  originated  a  small  distance 
aft  of  its  originally  calculated  position.  That  this  should  apply  to  the 
fluctuating  part  of  the  resistance  is  obvious  from  equation  (2.20),  which 
shows  that  the  smoothly  varying  part  of  the  resistance  is  a  function  of 
depth,  but  not  of  horlsontal  position.  The  extension  of  lnul9s  result  to 
the  resistance  equations  discussed  above  would  require  that  for  purposes  of 
calculating  the  interference  terms  the  positions  of  the  negative  sources 
(sinks)  at  the  after  end  of  the  hull  would  be  shifted  further  aft  by  a  small 
amount.  This  would  affect  the  coefficient  q^  in  equation  (2.22).  While 
this  has  not  been  demonstrated  by  experiment,  it  appears  to  be  a  reasonable 
extrapolation  from  Inui*s  findings. 

In  addition,  to  correct  the  scale  of  the  resistance  (slightly  in  error 
because  of  the  finite  height  of  the  waves) ,  Inui  multiplied  all  wave  heights 
as  computed  by  a  correction  factor  y<l.  This  would  translate  in  the  case 
of  our  equations  to  a  multiplication  of  each  mr  and  mg  by  the  factor  y  .  He 
found  further  that  the  wave  height  of  the  stern  waves  was  reduced  by  an 
additional  factor  |3 1 ,  and  this  would  require  that  all  sources  and  sinks  at 
the  after  end  of  the  ship  be  multiplied  by  a  factor  ft'  <  1.  Finally,  for 
purposes  of  calculating  the  interference  terms  only  he  reduced  the  amplitude, 
of  the  bow  wave  by  another  factor  a*  <  1 .  This  would  be  equivalent  to  multi¬ 
plying  the  amplitudes  of  sources  and  sinks  at  the  forward  end  of  the  ship  by 
the  factor  a?  in  the  calculation  of  the  Interference  terms  by  equation  (2.22). 
This  is  to  account  for  the  fact  that  the  resistance  equations  were  derived 
on  the  assumption  that  the  waves  could  propagate  over  the  entire  surface  of 
the  water,  but  they  are  actually  prevented  from  moving  aft  from  the  bow  through 
the  water  occupied  by  the  ship.  Therefore  some  part  of  the  bow  wave  is  less 
efficient  than  theory  would  predict  in  interfering  with  the  stern  wave. 
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While  the  theoretical  derivation  of  the  four  parameters  discussed  here 
is  sketchy!  they  have  worked  when  applied  to  bow  wave  and  stern  wave  ampli¬ 
tudes.  Since  any  ship  can  be  described  by  sources  placed  toward  the  bow  and 
sinks  placed  toward  the  stern!  it  appears  reasonable  to  try  to  transfer  the 
results  obtained  with  wave  amplitudes  to  the  formulation  using  sources. 
Certainly  the  case  of  a  hull  described  by  a  single  source  placed  forward  and 
a  single  sink  placed  aft  should  provide  a  satisfactory  candidate  for  the 
transfer:  the  bow  wave  starts  at  the  source ,  and  is  proportional  to  its 
intensity;  the  stern  wave  starts  at  the  sink,  and  is  also  proportional  to 
its  intensity.  Application  of  the  correction  factors  Inul  derived  for  the 
bow  and  stern  wave  amplitudes  would  therefore  apply  directly  to  the  strength 
of  the  source  and  sink  respectively!  and  the  correction  for  the  location  of 
the  stern  wave  would  apply  to  the  location  of  the  sink,  and  the  correspondence 
would  be  one-to-one. 

Since  the  calculated  curve  of  wave  resistance  vs.  speed  can  be  made 
identical  with  the  observed  curve  over  the  entire  range  of  speed  by  the 
introduction  of  only  four  parameters  (three  if  we  consider  that  only 
provides  a  proportionality  correction  between  the  curves  and  does  not  affect 
the  shape)  9  it  seems  reasonable  to  use  these  parameters  in  investigating  the 
results  of  small  changes  in  the  shape  of  a  known  ship.  That  is ,  if  we  can 
correct  the  theoretical  curve  for  a  known  set  of  sources  and  sinks  by  the 
use  of  these  parameters  so  that  it  is  identical  with  the  experimental  curve , 
then  it  is  probable  that  the  change  in  wave  resistance  we  calculate  for  small 
changes  in  the  magnitudes  of  these  sources  and  sinks  will  be  correct. 
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5.  RESISTANCE  OF  A  SEMI -INFINITE  STRUT  DESCRIBED  BY  A 


DISTRIBUTION  OF  DIPOLES  ON  THE  CENTERLINE  PLANE 
Om  form  which  hat  been  given  a  great  deal  of  investigation  is  the  smni- 
infinite  strut  which  can  be  described  by  a  distribution  of  dipoles  over  the 
centerline  plane  of  the  hull.  This  form  has  a  particularly  simple  wave* 
making  resistance  formula,  and  therefore  is  attractive  for  investigation. 

In  addition,  a  plot  of  the  density  of  dipole  moment  along  the  axis  of  the 
ship  looks  like  a  plan  view  of  one  side  of  a  ship  and  so  permits  easy  visualisa¬ 
tion  of  the  meaning  of  a  particular  distribution.  The  dipole  distribution  is 
usually  a  reasonable  equivalent  to  the  hull  shape.  However,  there  are  excep¬ 
tions  ;  for  example ,  Karp ,  Kotik  and.Lurye  found  a  case  where  the  dipole  density 
became  infinite  at  the  end  of  a  hull,  but  the  actual  width  of  the  hull  remained 
finite  and  the  shape  smooth  £ 6 j.  ' 

If  we  start  with  equations  (2.14)  and  (2.15)  and  substitute  sec20  -  cosh2u, 
we  get  the  following  expressions  for  the  case  k  ■  0  with  discrete  sources  mf: 

R  -  16ff^K02  T  (I2  +  J2)  cosh2u  du  ...(2.25) 


I  -  21  mr  cos(K0hr  cosh  u)  e"Kofrcoal,2u 
r-1 

J  -  £  m  sin(K0h  cosh  u)  e"Kofrco*h2u 
r-1 


...(2.26) 


We  may  now  take  any  source  and  an  equal  sink  (l.e.  a  source  of  negative  sign), 
the  sink  aft  of  the  source  by  a  distance  Ah,  and  combine  them  to  provide  a 
dipole  moment  mr  Ah.  The  value  of  this  product  we  can  call  ,M.  If  the  hull 
is  described  by  a  large  number  of  such  sources  and  sinks ,  it  is  possible  to 
combine  all  the  sources  and  sinks  at  any  one  depth  into  a  set  of  dipoles  with 
axes  oriented  in  the  direction  of  motion,  that  is  in  the  x-direction. 
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To  make  this  explicit,  let  us  assume  that  ve  have  a  set  of  n  discrete 
sources  at  a  depth  £  beneath  the  surface,  positioned  along  the  centerline 
plane  of  the  ship.  Since  the  x^dlrection  is  oriented  so  that  it  is  positive  > 
in  the  direction  of  motion  of  the  ship,  ve  will  assign  the  index  r  so  that 
the  value  r-1  corresponds  to  the  most  forward  source,  that  is,  the  one  with 
the  largest  value  of  h  for  its  x-coordinate .  The  index  will  increase  by 
one  unit  as  ve  come  to  each  source  aft  of  this  one.  Since  the  set  of  sources 
must  describe  a  closed  body,  it  is  necessary  that  the  sum  of  their  strengths 
be  zero* 

n 

Cmr  -  0  .  ,..(2.27) 

r-1 

We  will  define  the  distance  between  the  r***1  source  and  the  r+l8t  source  so 
that  this  distance  will  always  be  positive: 

m  hr  -  hr+1  ,..(2.28) 

It  is  now  possible  to  arrange  the  n  sources  r-l,...,n,  into  n-1  dipoles 
Mr  -  M(hp) ,  r  -  1 , . • . ,n-l . 

■  H(h^*)  -  m^4  h^  ,  h^ 1  -  h^  •  dh^/2 

Mj  m  110*2')  “  ^2  »  ^2*  *  ^2  * 

M3  -M(h3')  -(n^-H^-H*^) A h3  ,  h3*  «  h3  -ZJh3/2  ...(2,29) 


Mn-1  -  H(h,n.1)-(ml+...-Hnn.l)  4  \ml  ,  h^  -  h^- dh^/2 


The  tern  MQib^  terminates  the  sequence,  since  the  next  term  would  contain  the 
sum  (2.27),  equal  to  sero,  as  a  factor.  The  coordinate  h^9  is  the  x*coordlnate 
of  the  r**1  dipole  of  strength  H(hr*).  The  arrangement  Just  made  depends  on 
taking  the  source  m^  as  the  forward  member  of  the  first  dipole,  and  then 
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rap  lacing  tha  source  *2  by  the  sun  |^-n^+  (a^  +  .  The  portion  -Wj 

becoaee  the  after  aenber  of  the  first  dipole,  and  the  sua  (a^  +  ij)  becoaes 
the  forward  aenber  of  the  second  dipole.  Observe  that  the  dipole  M  is  , 
sctually  determined  by  the  strength  of  the  last  source  (which  will  have  negative 
strength),  aQ.  This  is  clear  since 

(Bj  +  b2  +  ...  -HsB.1)  -  -  b„ 

as  a  consequence  of  the  require  - 

sent  (2.27)  that  the  set  of  sources  describe  a  closed  body.  Fro®  this  It 
also  follows  that  the  set  of  dipoles  Is  the  same  whichever  end  of  the 

distribution  one  starts  with  In  making  up  the  source-sink  pairs  which  define 
the  dipoles. 

We  may  now  Investigate  the  special  case  where  Instead  of  discrete  sources 
n0  we  have  source  density  6  along  the  x-axis ,  and  the  function  6(h)  Is  Integrable. 
The  variable  h  Is  the  x-coordlnate  of  a  point  on  the  centerline  plane  of  the 
hull.  We  will  let  the  number  n  In  equation  (2  29)  become  Infinite  In  such  a 
fashion  that  maxdhr^>0.  Then  we  can  define  a  new  quantity ,  /i(h) ,  which 
we  will  call  the  dipole  density. 

/i(h)  -  11m  JM£h* )  .  ...(2.30) 

<jh-*0  dh 

In  the  limit  of  large  n  the  quantity  h*  will  coincide  with  the  quantity  h.  . 

Then  we  obtain  the  following  formula: 

h 

^i(h)  -  f  fdh  .  ...(2.31) 

The  lower  limit  of  the  Integral  In  (2.31)  Is  the  forward  end  of  the  source 
distribution. 

In  order  to  find  the  wave  resistance  of  such  a  distribution  we  must  return 
to  the  definitions  (2.25)  through  (2.29).  These  give  us  the  following  expres¬ 
sions  In  place  of  equations  (2.26): 
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I  ■  C  Mr  |coi^K0(h^<4hj/2)  co»h  uj-  cos|V0(h^  ^h^/2)  coah  9)  e«K0£rcosh*tt 


n-1 


■  I  gL.iln(K0h^  coah  u)  •inT^Ah^/i)  coah  ul  •~Kof*coih  u  ^ 

'■’am;  *•  J 

hm1  2 

-K0co*hu2I  Mr  sin(K0h^  cosh  u)  0”*®^°°'**  u  ,  for  &  mall.  6 hy  ■  4hr. 
r-1 

J  -  51  My  |  •in|lC0(h^+4 h^/2)  co»h  uj  -  sin|Ko(h^  -Ah^/2)  coah  ujje  *o^rco*^  u  m 

•  51  2Mr  cos  (KQhp  coah  u)  aln  f(^0  A  h^/2)coah  ul  e"K°^rC0*^  u 
r«l  XiTi 


A/ 

A/ 


n-l  _  f  .  2 

K0  cosh  u  Z  Mr  coa(K0hy  coah  u)  e”*°  *  u,  ^h£  small. 
r»l 


...(2.32) 


We  now  assume  that  there  are  many  layers  of  dipoles  of  identical  distribution 
but  different  depths,  so  that  we  can  separate  the  index  of  the  variable  f , 
the  vertical  position  on  the  centerline  plane,  from  the  variable  h,  the  fore- 
and-aft  position  on  the  centerline  plane.  Then  we  can  rewrite  equation  (2.32) 
as  a  double  sum.  We  assume  that  the  depth  of  the  hull  is  Infinite. 


n-l 

1  ~  -K  cosh  u  Z 

a-l 

n-l 

J  coah  u5Z 

a-l 


Z  Mg  sin(K0hJ  coah  u)  e‘K®frcosh  u 


r-1 

OP 


Z  M#  cos  (K0h|  cosh  u)  e~^°^rCOS^  u 


...(2.33) 


We  now  replace  the  Mg  with  /i(h,f)dh  df  -  /i(h)  dh  df  and  replace  tha  sums  with 
integrals.  Then  we  get  the  following  expressions  for  I  and  J: 
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I 


— Kq  cosh  u 


f  j  j»  (h)  «in(K0h  coih  u)  e"*®*  co*h2“ 
L  o 


dh  df 


1 


■  -(cosh  u)'1  )>»  (h)  sln(K0h  cosh  u)  dh 
L 

J  -  K  cosh  u  j  (h)  cos(K0h  cosh  u)  e'K°f  co*h2u  dh  df 
L  o 

-  (cosh  u)"1  J  fi  (h)  cos(KQh  cosh  u)  dh 
L 


>  ...(2.34) 


The  prime  hat  been  dropped  from  the  variable  h  because  la  patting  to  the  limit 
of  infinitely  small  4h^,  h  and  h1  become  coincident*  The  range  L  for  the 
integral  it  over  the  distribution  in  length  from  the  extreme  forward  source 
to  the  extreme  after  sink.  In  the  limit  of  continuous  dipole  moment  distribu¬ 
tion  the  range  of  source  positions  and  dipole  moment  positions  is  the  same. 

It  is  now  possible  to  calculate  the  resistance  of  this  dipole  distribution. 
We  square  the  two  expressions. 


I*  •  (cosh  u)“2  j  j  ji(h)  ji(h')  sin  (KQh  cosh  u)  sin(K0h*  cosh  u)  dh  dh'J 
L  L 

J2  •  (cosh  u)~2  J  ^  /i(h')  cos  (K  h  cosh  u)  cos(K  h'  cosh  u)  dh  dh' 

LL 


...(2.35) 


If  we  now  add  and  simplify  the  trigonometric  functions  we  get  the  following 
result: 


I2+J2  -  (cosh  u)"2  J  ^p(h)  ^(h1)  cosj^Oi-h1)  cosh  uj  dh  dh1  ...(2.36) 

L  L 

This  result  may  now  be  substituted  in  equation  (2.25)  to  obtain  the  resistance* 

K  -  lfty>Kc2  T  £  5>(h)  /»(h')  cot  T K0(h-h')  coah  uj  dh  dh'  du  ...(2.37) 
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We  aay  eliminate  ora  integration  by  uaa  of  the  identity 

5coajK0(h-h*)co3h  u]du  -  -(JT/2)  Y0(  |K0(h-h')|  )  -  -(y/2)  T0(K0|h-h'|  ) 

...(2.38) 

Than  we  get 

R  -  -  8fT2  p  K02  j  5y(h)  J»(h'>  VKolh"h'l  )  dh  dh,,  . .  .(2.39) 

L  L 

The  function  YQ(x)  is  Bessel's  function  of  the  second  kind.  Although  the 
function  Yq  has  a  singularity  at  the  origin,  its  Integral  is  a  bounded  function* 

For  positive  argument  Yq  Is  bounded  and  oscillates  from  positive  to  negative, 
gradually  diminishing  in  amplitude  with  Increasing  argument*  The  dipole  density 
p(h)  must  be  lntegrable  but  not  necessarily  continuous* 

6  a  RELATIONSHIP  BETWEEN  HULL  SHAPE  AND  SOURCE  DISTRIBUTION 
Since  the  resistance  equations  discussed  In  the  preceding  paragraphs  are 
all  based  on  calculating  the  resistance  of  a  collection  of  sources  and  sinks 
whose  closed  streamlines  outline  a  hull,  it  is  necessary  to  find  how  to  relate 
the  hull  shape  to  the  source  distribution.  In  principle  it  is  a  simple  matter 
to  find  the  hull  shape  corresponding  to  a  given  source  distribution,  but  it  is 
significantly  more  difficult  to  find  the  source  distribution  for  a  given  hull  shape 
One  of  the  assumptions  on  which  the  resistance  equations  are  based  is  that 
the  wave  height  is  small.  This  implies  that  the  flow  about  the  hull  is  not 
significantly  affected  by  the  existence  of  the  waves.  In  such  a  case  the  hull 
shape  can  be  described  by  the  closed  streamlines  on  either  side  of  the  plane 
of  symmetry  between  a  set  of  sources  and  their  mirror  image  in  the  free  surface , 
the  double  hull  so  described  moving  through  the  water  far  below  the  free  surface. 
Inul  £  3j  found  that  such  an  approximation  worked  for  his  models  up  to  a  Froude 
number  of  about  0.7.  Wc  can  therefore  write  the  velocity  potential: 


...(2. AO) 


n 

9  “  £  “a  <  i-  +  i_  > 

8ml  rls  r2s 

ru  ■  ^x-h^2  +  ^y_k^2  +  <z+f)2! 

r2»  "  £(*"h)2  +  (y-k)2  +  (z”f)2J 

It  is  clear  that  this  velocity  potential  will  provide  symmetry  in  the 

plane  s*0,  so  this  plane  may  be  considered  to  be  a  dividing  boundary  across 

which  there  is  no  flow.  This  also  follows  because  301  «0.  To  find  the 

9z  |z“0 

hull  form  corresponding  to  a  given  source  distribution  we  need  only  calculate 
0  from  equations  (2.40)  and  (2.41)  and  then  calculate  the  velocity  components 
from  0. 

u  *  -90  ;  v  •  -30  ;  w  *  -30  . 
fx  3y 

We  may  now  find  the  streamlines  by  a  numerical  integration  of  the  following 
differential  equation: 

dx  -  dy  «  dz  ,  ...(2.42) 

U  V  w 

Since  the  coordinates  move  with  the  ship,  this  equation  will  give  the  shape 
of  the  hull. 

Source  distribution  for  a  given  hull  form 
The  calculation  of  the  source  distribution  corresponding  to  a  given  hull 
form  is  a  much  more  difficult  matter.  John  L.  Hess  and  A.M.O.  Smith  have 
worked  out  a  way  of  doing  this  if  the  sources  are  distributed  on  the  surface 
of  the  hull  [9^.  With  such  a  technique  it  should  be  simple  to  provide  a  good 
estimate  of  the  wave-making  resistance  with  the  formulas  worked  out  earlier 
in  this  paper.  Unfortunately,  an  attempt  to  vary  the  strength  of  a  set  of 
sources  on  the  surface  of  the  hull  to  find  a  hull  of  improved  wave-making 
resistance  will  probably  run  into  difficulty  if  any  decreases  in  hull  volume 
are  permitted.  This  is  because  it  will  not  be  possible  to  define  a  hull 


surface  runs  inside  of  that  defined  by  these  sources  without  decreasing  their 
strength  to  zero.  Otherwise  we  will  calculate  a  hull  with  sources  outside 
its  surface  —  a  possible  arrangement  of  sources ,  certainly ,  but  one  which  , 
will  generate  Isolated  appendages  outside  the  main  hull* 

Notwithstanding  these  difficulties,  the  method  of  Hess  and  Smith  seems 
Important  enough  to  require  its  description  here.  We  may  assume  that  the 
given  hull  fora,  described  as  a  double  hull  in  an  infinite  fluid,  is  known 
and  given  by  the  following  equation: 

F(x,y,z)  ■  0.  •••(2.43) 

We  can  further  assume  that  the  velocity  of  the  ship  is  c  and  that  we  can 

write  the  velocity  potential  as  the  sum  of  a  term  cx  and  a  term  0*  which 

2  2  2 

vanishes  at  x  +  y  +  z  — *  «>*>  . 

0*0'+  cx  •  ...(2.44) 

We  may  also  describe  the  unit  normal  vector  directed  outward  from  the  hull 
by  if. 


...(2.45) 


How  the  normal  velocity  to  the  hull  surface,  deacribed  in  coordinates  which 
move  with  the  ship,  must  be  zero  at  the  hull  surface. 


-ag 

dn 


•  n 


grad  0 


0. 


...(2.46) 


We  may  substitute  the  relation  (2.44)  in  (2.46). 

n*  grad(0'+cx) 


[ 


-  ^  (0*  +  cx) 

5n 


12 

dn 


S 


L  ••[ 

£n*grad  0'J 


J 


-  0 


F-0 


*  -  n*i  c 


...(2.47) 


We  may  now  assume  that  the  surface  of  the  double  hull  is  covered  by  a  surface 
source  density  distribution  6.  Then  we  can  write  the  velocity  potential  as 
fallows: 
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...(2.48) 


0'(x,y,O  -  <£$  6(g)  dS 

S  TCPTqT 

where  r(P,q)  Is  the  distance  from  the  Integration  point  q  on  the  surface  to 
the  field  point  P  with  coordinates  (x,y,c)  where  the  potential  0*  is  being 
evaluated.  Kellogg  fioj  has  shown  that  the  normal  derivative  of  09  at  the 
point  p  on  the  surface  of  the  distribution  can  be  written  in  the  following  way: 

1  1  6(q)  dS  ...(2.49) 

*(P »q) J 

We  can  substitute  this  into  (2.47)  to  obtain  an  integral  equation  for  the 
distribution  6. 

2tt6(p)  -  55  iL  [  l 

S  dn  [r(p,q) 

Observe  that  ^i(p)  is  the  unit  normal  vector  at  F  *  0  and  the  variable  r(p,q) 
is  the  straight  line  distance  between  two  points  p  and  q  on  the  surface  of 
the  hull. 

Equation  (2.50)  has  been  Integrated  numerically  to  provide  not  a  source 
density  but  the  source  strengths  at  a  finite  number  of  points  on  a  double 
hull.  This  set  of  sources  provides  a  close  approximation  to  the  flow  about 
that  hull.  With  this  result  it  is  possible  to  calculate  the  wave-making 
resistance  directly  from  equations  (2.18)  to  (2.22)  -•  using,  of  course,  only 
those  sources  which  are  on  the  half  of  the  double  hull  which  is  submerged 
when  it  operates  as  a  surface  ship. 


6(q)  dS  *  n*  i  c 


...(2.50) 


-  -2 n  6(p) 


dn 
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111.  Methods  Which  Have  Been  Used  To  Find  Foma  Of  Decreased  Wave 
Resistance 

It  is  possible  to  draw  a  number  of  conclusions  from  the  form  of  the  > 
resistance  equations  and  from  the  Integrals  which  describe  forms  produced 
by  a  ship.  Beyond  this ,  Karp ,  Kotik,  and  Lurye  Q6 J  found  the  solution  to 
the  problem  of  the  form  of  that  semi-infinite  strut  of  minimum  wave-making 
resistance  which ,  between  the  ends  of  its  distribution  of  sources  and  sinks, 
can  be  described  by  an  integrable  distribution  of  dipole  moment.  The  first 
step  toward  a  method  of  steep  descent  was  made  by  Hogner  [ll]  in  1936,  but 
never  carried  through  to  completion.  Several  successful  attempts  to  find 
the  best  form  among  a  family  which  can  be  described  by  certain  polynomials 
subjected  to  various  constraints  have  been  made  —  all  of  these  using  some 
variation  of  the  Rits  method,  and  all  based  on  some  integrals  calculated  by 
Welnblum  £5j  ,  [12}  ,  fl3] .  A  more  striking  result  has  been  obtained  by 
Takao  Inul  [4]  by  determining  the  form  of  the  waves  produced  by  a  ship  and 
then  adding  to  the  hull  appendages  which  will  produce  roughly  similar  waves 
of  opposite  phase.  Finally,  certain  conclusions  have  been  drawn  by  Inui  as 
to  the  effect  of  discontinuities  in  the  source  distribution  or  any  of  its 
derivatives  on  the  waves  produced  by  a  serai-infinite  strut.  Some  further 
conclusions  can  be  drawn  from  these  as  to  what  kind  of  forms  cannot  have  very 
small  wave-making  resistance  without  being  very  small  themselves.  All  of  these 
matters  will  be  treated  in  turn  in  this  section. 

1 .  CONCLUSIONS  WHICH  CAN  BE  DRAWN  FROM  THE  FORM  OF  THE  RESISTANCE 
EQUATIONS  AND  OF  THE  INTEGRALS  DESCRIBING  THE  WAVE  FORMS 
Effect  of  Depth  on  Resistance 

One  of  the  more  obvious  conclusions  which  may  be  drawn  from  the  form  of 
the  equations  is  that  if  the  sources  describing  a  hull  can  be  submerged  deeply 
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enough  the  wave-making  resistance  can  be  made  as  small  as  we  please*  This 
follows  from  equation  (2.13)  which  contains  the  term  exp(-Kofsec2©)  ,  where 
K0-g/c2  and  f  is  the  depth.  On  the  other  hand,  the  resistance  is  certainly  ' 
non-negative,  since  equation  (2.12)  Is  the  integral  of  a  sum  of  squares 
times  a  positive  definite  trigonometric  function* 

A  further  conclusion  of  the  same  sort  can  be  drawn  by  examining  equations 
(2*19)  and  (2*20)*  These  equations  show  that  the  smoothly  varying  portion  of 
the  wave-making  resistance  is  small  for  very  low  K0#  Increases  to  a  maximum, 
and  then  decreases  again*  For  a  proof  see  appendix  1*  Since  K0"g/c2,  it  is 
clear  that  the  largest  value  of  the  resistance  occurs  at  lower  speeds  for 
small  f  (small  depth)  and  at  higher  speeds  for  large  f (large  depth).  This 
means  that  for  high  enough  speeds  the  contribution  to  the  smoothly  varying 
part  R^  of  the  wave-making  resistance  from  parts  of  the  hull  corresponding 
to  sources  very  close  to  the  surface  will  approach  zero.  The  interaction 
term  R^  displays  a  similar  behavior,  although  it  is  clear  that  the  inter¬ 
action  between  a  shallow  source  and  a  deep  one  will  reach  its  maximum  absolute 
value  at  a  speed  which  is  intermediate  between  the  speed  for  the  maximum  of 
the  smoothly  varying  term  for  the  shallow  source  and  that  for  the  deep  source. 
The  general  conclusion  which  can  be  drawn  from  these  facts  is  that  the  part 
of  the  hull  near  the  waterline  can  be  designed  for  relatively  low  resistance 
at  low  speeds,  and  the  part  of  the  hull  deep  in  the  water  for  relatively  low 
resistance  at  high  speeds,  and  the  interaction  terms  between  shallow  and  deep 
sources  can  be  used  to  improve  the  behavior  at  high  speeds* 

Effect  of  Symmetry  Fore  and  Aft 

Another  conclusion,  this  on  old  and  frequently  misstated  one,  is  that 
which  relates  to  the  effects  of  fore-and-aft  syrmetry.  Suppose  that  wo  have 
a  dipole  moment  density  distribution  such  that  the  after  half  of  the 


29 


distribution  is  a  reflection  of  the  forward  half  (i*e.t  if  the  origin  of  x 
is  at  the  midships  section,  p(x)  is  even)*  Then  the  function  1  in  equation 
(2*34)  Is  zero  and  the  resistance  of  (2*25)  reduces  to  ^ 

R  -  l6»pK02  }  J2  coah2u  du  ...(3.1) 

o 

How  if  we  add  to  /i(x),  which  we  have  assumed  to  be  even,  any  odd  function  of 
x  there  will  be  no  change  in  J  but  1  will  become  different  from  sero*  Then 
ve  will  have  to  add  to  J*  in  (3*1)  a  positive  definite  quantity  and  so 
will  Increase  the  resistance*  Since  the  odd  function  added  to  /u(x)  will  not 
change  the  total  dipole  moment  (which  is  roughly  proportional  to  volume) 9  it 
follows  that  the  hull  described  by  an  even  function  u(x)  has  less  resistance 
than  the  hull  of  roughly  the  same  volume  described  by  that  even  function  plus 
an  odd  function  of  x*  This  is  not  the  same  thing  as  a  statement  that  a  hull 
symmetrical  fore-and-aft  has  less  resistance  than  one  which  is  asymmetrical  -- 
it  is  extremely  easy  to  devise  a  symmetrical  hull  with  more  resistance  than 
any  given  asymmetrical  hull  of  the  same  volume*  The  demonstration  of  this 
paragraph  can  be  generalized  to  hulls  described  by  an  arbitrary  symmetrical 
fore-and-aft  distribution  of  dipoles  plus  an  arbitrary  antisymmetrical 
distribution.  We  need  only  use  equation  (2*32)  Instead  of  equation  (2*34). 

It  is  possible  to  argue  that  the  experience  of  practical  ship  designers 
does  not  bear  out  the  conclusion  that  the  hull  of  least  resistance  is 
symmetrical  fore-and-aft*  This  can  also  be  demonstrated  by  appealing  to  the 
semi -empirical  parameters  used  by  Inul  to  correct  his  theoretical  resistance 
curves  to  coincide  with  those  found  from  experiment.  For  example,  he  found 
that  he  had  to  reduce  the  effect  of  the  sources  at  the  extreme  after  end  of 
the  ship  by  a  factor  p* <  1  in  order  to  bring  the  calculations  into  accord 
with  observation*  This  is  because  the  after  end  of  the  ship  is  relatively 
less  effective  in  making  waves  than  the  bow.  Wc  can  do  this  in  our  formulation 
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by  multiplying  each  source  at  the  after  end  of  the  ship  by  p1  and  then  repeat¬ 
ing  the  derivation  of  equation  (2.32).  Then  the  source  distribution  which 

-  > 
makes  r  »  0  will  be  symmetrical  fore  and  aft  in  the  moments  derived  from 

sources  multiplied  (at  the  after  end)  by  p1;  and  this  will  certainly  not  be 

symmetrical  fore-and-aft  in  the  actual  dipole  moments,  Mr» 

Effect  of  Length 

It  has  long  been  known  empirically  that  the  resistance  per  unit  of 
volume  tends  to  decrease  as  the  ship  is  made  longer.  This  is  not  true  for 
all  small  changes  in  length:  for  example,  if  the  length  is  such  that  the 
transverse  bow  and  stern  waves  tend  to  cancel  each  other,  increasing  the 
length  a  small  amount  may  cause  them  to  reinforce  and  increase  the  resistance. 
On  the  other  hand,  for  large  changes  in  length  the  general  tendency  will  be 
determined  by  the  smoothly  varying  portion  of  the  resistance  expression. 

But  we  showed  in  equation  (2.20)  that  the  smoothly  varying  portion  of  the 
resistance  is  a  function  only  of  the  strength  of  the  sources  and  of  their 
depth,  not  of  their  horizontal  position.  On  the  other  hand,  if  we  have  a 
single  source  forward  of  an  equal  sink  moving  through  the  water  at  a  constant 
speed,  the  set  of  closed  streamlines  which  they  generate  will  increase  in 
volume  monotonlcally  as  they  are  moved  farther  and  farther  apart.  It  follows 
from  this  that  the  longer  is  a  hull  described  by  a  single  source  and  a  single 
sink  of  constant  intensity,  the  smaller  will  be  the  resistance  per  unit  volume 
if  we  consider  only  the  part  of  the  resistance  which  varies  smoothly  with 
speed*  The  same  result  can  be  obtained  if  any  assembly  of  sources  and  sinks 
which  produces  a  closed  set  of  streamlines  is  moved  farther  apart  in  the 
x-direction. 
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Effect  of  Breadth 


For  a  very  thin  hull  it  has  been  accepted  that  the  source  density  can 

* 

be  approximated  by  equation  (2.24): 

6  -  -  c_  ik  . . .(2.24) 

2w  dh 

It  follows  that  if  we  multiply  the  source  density  6  by  the  parameter  a  we 
multiply  the  breadth  of  the  hull  by  a: 

h 

k  -  -  2?ra  C  6  dh  ...(3.2) 

c  bow 

The  minus  sign  results  from  the  fact  that  the  axis  of  x  is  positive  forward. 

Now  if  we  substitute  a6*  for  6  in  equation  (2.15) ,  we  get 
P(a6*)  -  aP(6f) 

R(a6’)  *<  l6rpKo2  $  Jp2(a6')  +  Q2(a6’)J  sec3©  d9  - 
-  l6iyaK02  a2  jP2(6’)  +  Q2(6')J  sec3©  d©  -  a2R(6').  ...(3.3) 

From  (3.2)  and  (3.3)  it  is  clear  that  the  resistance  of  very  thin  hulls  is 
proportional  to  the  square  of  the  breadth.  Inui  £  3j  and  Hess  and  Smith  £9^| 
have  shown,  however,  that  for  hulls  of  finite  breadth  this  relationship  is 
Incorrect.  No  general  relationship  has  been  worked  out  for  such  hulls,  and 
it  is  clear  from  the  plots  that  have  been  obtained  that  any  such  relationship 
would  be  very  complicated. 

Effect  of  Bow  Shape  on  Resistance 

If  the  source  distribution  describing  the  hull  consists  of  a  distribution 
over  the  centerline  plane  of  the  ship  in  such  a  fashion  that  it  is  not  a 
function  of  depth  but  is  a  continuous  function  of  length  and  has  continuous 
derivatives  of  all  orders,  it  is  possible  to  calculate  the  effect  of  the 
shape  of  the  bow  on  resistance.  Let  us  consider  equation  (2.9a)  and  substitute 
for  the  assemblage  of  discrete  sources  a  source  distribution  6(h,0,-f)  *  6(h) 
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od  the  centerline  plane .  Then  we  may  Integrate  with  respect  to  f  and  take  f 


to  Infinity  (infinite  depth). 
irll  h—L 


y  4K0  J  J  5  i(l>)  co»|i0(x-h)  see  #j  coa  ^Kpysec*®  ilnoj 
x  exp(-K0fsec^8)  sec^O  d®  dh  df  - 


W/I  h--L  r  i  r  9  1 

4  J  f  6(h)  coslK0(x-h)  sec®  I  cos  |Koyaecz0  ain®|  sec®  d®  dh 

•  -m  h-o  L  J  L  J 


...(3.4) 

It  is  now  possible  to  integrate  by  parts  with  respect  to  h  so  that  we  obtain 
a  series  in  6(h)  and  its  successive  derivatives,  evaluated  at  h  ■  0  (the  bow) 
and  h  ■  -L  (the  stern). 

H  fL  6(h)  cos  Ik  (x-h)  sec®]  sec®  dh  -  -1  f 6(h)  slnflC  (x-h)  aecoll  + 

h-o  L  1  5o[  L  Jj  h-o 

h“-L  p  -j 

+1  J  6'(h)  sin/K  (x-h)  sec  ®|  dh  - 
Ke  h-0  L  J 

-  1  6(0)  ain(K0x  sec®)  -  l  €(-L)  sln^a+x)  sec®)  + 

K©  Ko  L 


+  .,-J— 

Ko‘sec0 


h“-L  hw-L 

-T 


h-0  Kp' 


sin(K  x  secO) 
o 


C08(Kqx  sec®) 


js 1  (h)  coajKp(x-h)sec®j^ 

-  1  f  6(0)  -  6m(0)  +  6lv(0)  -  . . . 1 

K©  j  K0^sec^®  K^sec*®  J 

-  1  f  6'(0)  -  6''*(0)  +  ....1 

Kp^KpSec®  V’sec*®  J 

-if  6(-L)  -  6"(-I.)  +  6tv(-L)  -  ...(  sln|K  (L+x)scc®l 

Kpl  K^sec2®  Kp^sec4®  V  L  J 

+  1  f  6*(-L)  -  6,,,(-L)  +  ....?  cos/k  (L+x)sec®l 

*o\Ko0CC®  Ko3Sec-*0  J  L  J 

»  S(0,G)  sin  (Kcx  sec  9)  +  C(0,9)  cos  (K0x  sccO) 

-fs(-L,®)  sinjK0(x+L)sec®|  +  C(-L,0)co3  ^Ko(x+L)ncc0J  C 

l*  33  ^ 


h--L  r  •» 

1 _  V  6"(h)  cos  K  (x-h)secOl  dh 

sec®  h-0  L  °  J 


...(3.5) 


where 


S(h, 

C(h 


®)  -1  { 


6(h)  -  £^hl  .  +  6lv(h) 


K^Zaec2©  K^tec^S 
.9)-  -1  /tf’Ch)  -  6*"  (h) 

*0^0' 


•] 


.(3.6) 


^  m--,  +  e  e  e  ^ 

|K0sec9  Veec^g  j 


We  can  now  write  equation  (3*4)  in  terms  of  S  and  C: 

r/r2  9  9 

t>«  4  S  8(0,8)  sln(K0x  eec20coa0)coc (K-y  eecz8sln8)  d8  + 

J  c  -Tf2 

H2  2  9 

+  4  5  0(0,9)  cos(K  x  eec40cos9)coa(K  y«ec20  alnO)  d6  - 

c  -T/2 

yjj2 

-4  5  S(-L,0)  ain|k0(x+L)  sec20eos0j  cos(Koysec20  alnO)  dO  - 

c  -Tfl2 

^j2 

-4  ^  C(-L,0)  cosTk  (x+L)sec29coscl  cos (Ky sec 2 8  alnO)  dO 

c  -Tf2  L°  Jo  ...(3.7) 

It  la  now  poaaible  to  write  the  wave-making  realatance  correapondlng  to  thla 

wave  forme  For  this  purpose  we  will  examine  only  the  first  two  terms  of  the 

expression  (3,7) ,  which  are  the  waves  generated  at  the  bow  of  the  ship*  If 

we  disregard  the  interference  between  bow  and  stern  wave,  the  resistance 

produced  by  the  bow  waves  can  be  written  as  follows  [3j: 
nfj  f  r  ^  r  , 

cos3 9  d9 


R  -  const 


S  |^s(o,e)p  +  [c(o,e)|2| 


...(3.8) 


For  very  low  velocities  KQ  is  very  large,  and  go  the  dominant  term  in  (3*8) 
is  the  first  term  in  the  expression  for  S(0,0)  in  (3.6)  —  that  is,  the 
source  density  at  the  bow.  Therefore  the  smaller  this  source  density  the 
smaller  the  resistance  for  very  low  speeds*  The  source  density  at  the  bow 
is  proportional  to  the  angle  of  entrance  at  all  points  except  at  the  extreme 
bow,  so  it  is  clear  that  for  very  small  speeds  the  wave-resistance  produced 
by  the  bow  waves  is  determined  by  the  entrance  angle  of  the  bow  -•  the  larger 
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the  angle,  the  larger  the  resistance.  However,  as  the  speed  increases  K0 
becomes  smaller  and  the  later  terms  in  the  expansion  for  S  and  C  become  import¬ 
ant,  and  when  KQ  is  of  the  order  of  unity  they  become  of  nearly  the  same  order 

of  importance  as  the  source  density  itself.  For  KQ  small  enough,  for  example, 

it  is  possible  to  have  the  first  and  second  terms  of  S(0,9)  become  of  equal 
magnitude  and  opposite  sign  so  that  they  eeneel  fei  values  of  0.  Then 

if  all  other  derivatives  (0,0),  k»l,3,4,...  are  zero  or  small  in  magni¬ 


tude  we  may  find  that  a  ship  of  small  resistance  will  have  a  large  source 
density  at  the  bow  and  a  small  density  farther  aft.  This  implies  a  bluff 
bow  with  large  entrance  angle. 

Non-Zero  Resistance  of  a  Continuous  Symmetric  Source  Distribution 
One  further  conclusion  may  be  drawn  from  equations  (3.4)  to  (3.8),  We 
may  take  advantage  of  the  fact  that  S(-L,0)  *  -S(0,Q)  and  C(-L,9)  *  C(O,0) 
for  a  hull  which  is  symmetrical  about  its  midship  section  (i.e.,  fore  and 
aft).  Then  we  have  the  following  expression  for  the  wave  height  ^  : 

•  4  /  5(0,8)  isin(K0x8CC28cos9)  +  ainfK0(x+L)sec2Qcos9l  •  cos(Koyscc203in8)  d0  + 

c  -7/2  (_  L  Jj 

+  4  J  C(0,9)|  coa(Koxaec20cos6)  -  cos|^K0(x+L)scc29cosel  V  cos  (K0ysec29sin8)  d9  = 

^  V 

-8  J  S(0,9)  sln|Ko(x+L/2)scc20cos0j  cos  jv Co(L/2)oec28cos9|  cos(Koyscc20sinQ)d8  ^ 

+  8  J  C(0,8)  ainjK0(x+L/2)aec2ecosGj  sin  j^Ko(L/2)sec20cos9j  co3(Koyscc20ainO)d0 

We  may  combine  the  factors  injK0(L/2)sec2ecosfl]  with  S(0,0)  in  the  first  term 


..(3.9) 


and  with  C(0,Q)  in  the  second  term.  We  can  write  then 
S^-L.0)  -  C(0,9)alnjKo(L/2)sec20co50j 
S2(-L,9)  -  S(0,9)cc.sjKo(L/2)scc20cos0j 


...(3.10) 
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It  follows  that  the  resistance  of  the  ship  is  given  by  the  following  expression: 


7r/2 

R  ■  const  1 

rs.(-L,o)+ 

s,(-L,ejj  2  co.3e  do 

...(3.11) 

o 

L 1  2 

*  2  J 

Since  cos  ©  is  real  and  non-negative  for  the  entire  range  of  integration , 

and  the  squared  sum  in  the  Integrand  is  also  real  and  non-negative ,  -S^ 

for  SpS2fOt  R  can  be  made  zero  only  if  and  $2^1,0)  sre  separately 

2  2 

zero  for  all  0  in  the  range  0  <  0  <  tH2.  But  this  can  be  true  only  if  6(0), 
6*(0),  6M(0),...  are  separately  xero,  since  the  series  in  (3.6)  cannot  other¬ 
wise  be  sero  for  all  0  in  this  range.  It  follows  that  the  wave-making 
resistance  of  an  infinitely  deep  strut  which  is  symmetrical  fore-and-aft, 
which  is  described  by  a  continuous  distribution  of  source  density  on  the 
centerline  plane,  which  is  not  a  function  of  depth,  and  which  has  continuous 
derivatives  of  all  orders  cannot  be  sero.  Further,  it  follows  from  a  previous 
result  that  if  a  continuous  fore-and-aft  anti -symmetrical  distribution  of 
dipole  moment  density  is  added  so  that  the  hull  loses  its  fore-and-aft 
syametry ,  the  resistance  can  only  be  Increased.  This  result  says  nothing 
concerning  the  resistance  of  an  assemblage  of  discrete  sources  and  sinks. 
Effect  of  a  Discontinuity  in  Source  Distribution  or  Any  Derivative 
Let  us  assume  in  the  derivation  of  equation  (3.5)  that  the  k**1  derivative 
has  a  discontinuity  at  h«h^.  This  will  give  rise  to  a  term  of  the  form 
6(k)(hk+)  -  6(k)0y)  , 

O 

so  that  in  addition  to  the  wave 

pattern  described  by  equation  (3.7)  we  have  a  wave  described  by  the  following 
equation: 


7f/2 

+  4  P  1  6^) (hk+)-6^) (hk-)  (sin  ■»  "1  o 

'  ~  c  -ir/2  K0 - ^  k  gke -  I  cos  K0(x4hk)sec2Gcose|  cos(K0yoecZOslnO)  dO 


...(3.12) 
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where  the  sine  or  cosine  is  taken  depending  on  whether  the  integer  k  is  even 
or  odd,  and  the  sign  is  positive  for  k  -  4n,  k  -  4n-l,  n  »  0,1,2,...,  end 
negative  otherwise.  This  explains  the  existence  of  the  so-called  shoulder  j 
waves  of  a  hull,  which  start  at  such  discontinuities  as  the  connection  between 
the  fair  bow  section  and  the  parallel  Middle  body.  Since  these  additional 
waves  occur  at  a  point  where  they  cannot  in  general  be  fully  canceled  by 
waves  originating  at  the  bow  and  stern,  they  will  add  to  the  wave  resistance. 

2.  APPLICATION  OF  THE  CALCULUS  OF  VARIATIONS  TO  FINDING 
HULL  FORMS  OF  MINIMUM  RESISTANCE 

An  apparently  reasonable  approach  to  the  problem  of  finding  a  ship  of 
mlnlmua  wave  resistance  is  to  apply  the  classical  calculus  of  variations. 
Unfortunately,  the  failure  to  distinguish  between  the  hull  form  and  the 
dipole  density  which  generates  it  together  with  the  general  difficulty  of 
the  problem  prevented  a  solution  being  obtained  for  many  years.  Sretenskij 
M  published  a  paper  in  1935  in  which  he  demonstrated  that  in  certain 
cases,  at  least,  there  was  no  solution.  Even  today  the  problem  has  only  been 
solved  for  the  case  of  a  strut  of  infinite  depth. 

Karp,  Kotik,  and  Lurye  £l5j  solved  the  problem  in  1959.  They  first  con¬ 
sidered  a  dipole  moment  distribution  uniform  in  depth  on  an  infinitely  deep, 
centerline  plane  of  length  L,  and  assumed  that  the  dipole  moment  density 
vanished  at  the  bow  and  stern  and  had  an  integrable  first  derivative  with 
respect  to  x.  They  then  were  able  to  prove  that  if  an  additional  requirement 
were  imposed  that  the  dipole  moment  per  unit  depth  be  held  constant,  the 
Integral  equation  resulting  from  the  application  of  the  classical  calculus 
of  variations  had  no  solution.  They  reasoned  that  the  cause  was  that  the 
class  of  functions  chosen  for  the  dipole  moment  density  was  too  restrictive , 
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and  then  tried  a  more  general  aet  of  function*.  This  choice  of  function* 
required  that  the  velocity  potential  not  become  infinite  at  finite  points 
of  the  (x9y)  plane  other  than  the  centerline  plane  on  which  the  dipoles  are 
distributed  and  that  the  integral  of  the  dipole  moment  over  the  length  L 
remain  constant  per  unit  depth.  With  this  extended  choice  of  functions  they 
were  able  to  find  a  solution  for  the  dipole  moment  distribution  which  pro¬ 
vided  a  form  of  minimum  wave-making  resistance. 

In  order  to  sketch  the  method  employed  by  Karp,  Kotik,  and  Lurye  their 
notation  will  be  converted  to  that  used  elsewhere  in  this  paper.  We  will 
start  by  assuming  that  the  velocity  potential,  as  before,  can  be  written  as 
the  sum  of  a  term  corresponding  to  the  flow  in  the  absence  of  the  ship  and 
a  second  term,  which  vanishes  at  infinity,  corresponding  to  the  disturbance 
caused  by  the  distribution  of  dipoles  along  the  centerline  plane  of  the  hull. 

The  axis  of  x  will  be  positive  lu  the  direction  of  ship  motion. 

0  »  0'  +  cx  • ..(2.44) 

The  dipole  distribution  will  be  taken  along  a  strip  of  infinite  depth  and 
length  L.  We  will  lose  no  generality  if  we  assume  that  the  length  L  of 

the  distribution  is  unity.  We  can  always  correct  to  other  lengths  by  Froude's 
law  of  similitude.  Since  we  have  already  demonstrated  that  if  a  form 
symmetrical  fore-and-aft  can  be  found  whose  resistance  is  minimised  any 
addition  of  an  anti -symmetrical  form  to  it  will  only  Increase  the  resistance, 
we  may  also  assume  that  the  hull  form  is  symmetrical  about  the  midships  section. 
We  will  write  the  dipole  moment  distribution  in  terms  of  a  new  variable 

^-h  +  l/2  ...(3.13) 

Then  we  will  choose  a  value  for  the  Integral  of  the  dipole  moment  over  the 
length  of  the  distribution.  If  this  is  held  constant,  then  the  volume  per 
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wait  depth,  which  is  approximately  proportional  to  the  dipole  moment,  will 
be  held  constant  to  first  order.  He  will  let  this  integral  be 


1/2  _  . 

S  1»<«  i\  ■  A >  0.  ...<3.14) 

-1/2  *  > 

The  quantity  is  the  dipole  moment  per  unit  depth  and  unit  length.  The 

other  restriction  to  be  placed  on /*(<£)  is  that  the  following  integral  for 
the  velocity  potential  be  well-behaved: 


1  /2  w  y 

•  const  )  ^  ^»<T)  *C(x,y,t;  S,0,-f)  d  ^  df  ...(3.15) 

°  -1/2  #1 


The  function  G(x,y  ,s;  ,0,-f)  is  the  velocity  potential  of  a  unit  source  and 

can  be  derived  from  the  first  term  of  equation  (2.4)  by  letting  m  *  1  and 
substituting  equation  (3.13). 

G(xty9s; J  ,0,-f)  « _ 1 _  ...(3,16) 

l  -  +  y^  +  (z+f )* 


The  partial  derivative  with  respect  to  ^  converts  G  into  the  Green's  function 
for  a  dipole  with  axis  oriented  along  the  x-axls.  The  restriction  placed  on 
^(£)  in  addition  to  (3*14)  is  that  0*  be  finite  everywhere  except  on  the 
plane  on  which  the  dipoles  are  distributed,  and  that  it  go  to  zero  as 
x2  +  y2  +  z2  — The  quantity  x  is  still  reckoned  from  the  forward  end 
of  the  dipole  distribution. 

In  part  11  of  this  paper  we  proved  that  the  resistance  of  an  assemblage 

of  dipoles  can  be  written  in  the  following  way  (see  equation  (2.39)): 

1/2  1/2 

R.-Sir2^2  $  5  VKon-11>  dl  dV  • 

-1/2  -1/2  *  ' 

He  may  regard  the  coefficient  outside  the  integral  as  a  constant  and  combine 
this  equation  with  (3.14)  to  write  the  variational  equation: 
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i ■ «  i  j  |V<5>  * •?<!>]  ^Yf\ >[,<;•> + * 

+  xj|f  [/•(j)  *  4J  *  AJ  ...(3.18) 

We  may  now  find  the  Euler  equation  by  taking  bj  and  letting  e  — >0. 


a 

d* 


«— >o 


h* 

1/2  1/2  .  , 

2c  j  y}’>  5  -j’  >  dldJ*  + 

-1/2  -1/2  ’ 

+1  *f  i(J’)  <*|  '  “  o 

-1/2 


...(3.18a) 


From  this  we  can  write 


...(3.18b) 


l/2  f  1/2  „  ,v  „  /  v  1  v 

^  W>  2C  J  P<J>  Yo^o  f -J'p f  “°  • 

-1/2  /  -1/2  J  } 

But  the  function  ^(^')  is  arbitrary,  so  it  follows  that  the  quantity  in 
braces  is  zero  for  all  This  is  the  Euler  equation  for  the  problem.  The 

linear  nature  of  the  equation  in  braces  makes  it  possible  to  solve  the  equation 
for  arbitrary  A  and  then  normalize  the  solution  by  multiplying  by  the  proper 
factor  to  satisfy  equation  (3.14),  which  we  regain  from  (3.18a)  by  taking 
bf  |  *  0.  Hence,  we  may  write  without  loss  of  generality, 

J«-^0 

1  1/2  ^  ....  i  .  . 

...(3.19) 


1/2 

$ 

■1/2 


and  solve  for  A1  ■  1  and  then  normalize  the  solution  to  (3.14).  Equation 
(3.19)  has  been  shown  [l5j  to  have  solutions  of  the  form 


)  -  K<n 


> 


joTtrj 


...(3.20) 
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where  g(  j£)  is  regular  in  the  range  -1/2  <  <  1/2,  and  not  aero  at  the 

endpoints.  In  consequence  the  density  of  the  dipole  moment  goes  to  infinity 
at  each  end  of  the  distribution*  It  turns  out,  however,  that  this  does  not  » 
result  in  the  velocity  potential  becoming  infinite*  We  may  show  this  by 
substituting  (3*20)  in  (3*15)*  Then  we  may  examine  the  portion  of  the 
velocity  potential  which  results  from  the  singularity  at  *  1/2.  We  will 
call  this  A  V  . 


r  111 

A  0'  « const  \  df 

°  1/2 


/2  v  * 
5  «<J>  H  ■—£ 

2-e  y  >/(l/4)  - 


l 


const  \  df  g(l/2-  *e)dG(l/2-«0)  \  df* 

-  const  \  g(l/2*c©)dC(l/2-*0,O>-f)  fsin"ll  -  sin-1(l-2i)/  df 

o  *5  L  -J 

But  »in"*(l)  -  sln’*(l-2«)  «  cos“*(l-2e)  ^  2^?  ;  and  so 


Af'^const  $  g(l/2-e0)2|G(l /2-e0,O ,-f ) (2>/«)df ,  O<0<1.  ...(3.21) 

• 

Equation  (3.21)  is  still  bounded  after  integration  with  respect  to  f,  and 
so  the  contribution  of  the  singularity  to  0'  is  also  bounded  and  the  condi¬ 
tion  on  ji(Jf  )  is  met. 

Karp,  Kotik 9  and  Lurye  have  solved  equation  (3.19)  numerically  for 
numerous  cases  and  have  plotted  the  results  as  streamlines  [6j.  The  closed 
streamlines  which  define  the  boundary  of  the  strut  extend  farther  forward 
and  farther  aft  than  the  source  distribution.  For  large  values  of  the  speed 
they  have  found  that  the  strut  of  minimum  resistance  has  a  cross-section 
shaped  rather  like  a  dumbbell  --  rounded  at  the  leading  edge,  then  narrow, 
then  expanded  and  rounded  at  the  trailing  edge.  For  very  large  speeds  (Froude 
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numbers  greater  than  about  0.65)  the  function  g  (^)  in  equation  (5.20) 
becomes  a  constant.  For  small  speeds  the  function  g(J)  *  minimum  at  the 
end  points  and  a  maximum  at  the  midships  section. 

There  is  a  simple  interpretation  which  will  make  the  infinities  at  the 
ends  of  the  dipole  distribution  appear  more  reasonable  than  they  might  seem 
at  first.  If,  for  example,  ve  started  at  the  forward  end  of  the  distribution 
with  an  Isolated  dipole  of  strength  H,  and  a  source  of  strength  m,  while  at 
the  after  end  of  the  distribution  we  had  another  isolated  dipole  of  strength 
M  and  a  sink  of  strength  -m,  we  would  have  a  set  of  singularities  which  are 
well  known  to  generate  a  reasonable  set  of  streamlines  in  two-dimensional 
flow.  If,  further,  we  identified  the  isolated  dipole  with  a  delta-function 
multiplied  by  a  strength  H,  then  we  would  have  a  condition  where  the  dipole 
density  became  infinite  within  the  bounds  of  the  delta- function,  but  the 
total  moment  remained  finite  and  equal  to  H.  The  distribution  obtained  by 
Karp,  Kotik,  and  Lurye  for  very  high  speeds  is  not  too  different  from  the 
condition  just  described. 


3.  USE  OF  THE  RITZ  METHOD  FOR  FINDING  HULL  FORMS  OF 
REDUCED  WAVE-MAKING  RESISTANCE 

The  bilinear  nature  of  the  resistance  integral  leads  in  a  natural  way 
to  the  use  of  the  Ritz  method  for  finding  ways  to  develop  improved  hull  forms 
It  has  been  used  successfully  for  this  purpose,  although  its  success  has  been 
limited  both  by  the  restricted  set  of  functions  employed  for  the  description 
of  the  ship  and  by  the  (to  date)  consistent  failure  to  distinguish  between 
the  distribution  of  dipole  moment  and  the  actual  shape  of  the  hull.  One  of 
the  more  striking  successes  was  the  design  of  a  hull  which,  when  towed  at  a 
Froude  number  of  0.5,  had  about  13  per  cent  less  wave-making  resistance  than 
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a  destroyer  of  similar  dimensions,  end  which,  even  after  the  other  components 
of  resistance  were  included,  had  about  8  per  cent  less  total  resistance  than 

> 

Che  destroyer  [i2j. 

In  order  to  see  the  problem  in  simple  terms,  we  will  start  with  equations 
(2.14)  and  (2*15).  Here  we.  have  a  source  density  6  which  generates  the  hull 
form.  If  ve  add  to  this  source  density  another  distribution  a6' ,  then  ve 
have  a  resistance  formula  which  is  a  quadratic  in  the  parameter  a. 

o  ^ 


n/2 

R  -  l6n(>K 02  J  £p(*)2  +  2aP(6)P(6')  +  a2P<«')2  +  Q(6)2  +  2aQ(6)Q(6')  +  0^(6*)^ sec3  9  < 


This  can  be  rewritten  as  follows: 

R  -  R(6 ,6)  +  2aR(6,6’)  +  a2R(6',6') 

where  the  functional  R(6,6')  is  defined  by 

0  r/2  r  ^  * 

R(M')  -  l6*pKo2  J  jp^Ptf1)  +  Q(6)Q (6')J  sec30  d9 


...(3.22) 


...(3.23) 


...(3.24) 


It  is  now  a  simple  matter  to  take  the  derivative  of  R  in  equation  (3.23)  with 

respect  to  the  parameter  a  and  find  the  optimum  amount  of  the  distribution  6’ 

to  add  to  the  given  distribution  6  by  setting  the  derivative  dR  »  0. 

da 

To  generalize  this  result,  we  need  only  follow  the  same  pattern  but  give 

each  of  arbitrarily  many  linearly  independent  source  distributions  # *  *^n 

a  coefficient  a, a  and  write 
I  i  n 


R  -zz 


n-l 


...(3.25) 


mZ~'a  R (6r,6r)  +  2  r  E  ara8R<*r»6*> 

r»l  9-r+l  r-1 

If  the  set  of  functions  ^6^  ,6^ , • • is  also  a  complete  set,  then  its  members 

may  be  uaed  to  describe  an  arbitrary  dipole  moment  distribution  on  the  interval 

on  which  they  are  defined.  It  is  a  simple  matter  also  to  add  constraints.  We 

merely  write  them  in  the  form 
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%*••)  *  ®  ...(3.26) 

*nd  then  apply  the  method  of  Lagrange  multipliers  to  find  a  stationery  value 
of  K  subject  to  these  constraints.  To  do  this  we  write  ' 

|  -  R+lXlf1  ,1-1,...  ...(3.27) 


where  R  is  defined  in  equation  (3.25),  and  then  take  the  partial  derivative 
of  J  with  respect  to  each  ay  and in  turn  and  set  it  equal  to  sero. 


The  functions  6r  which  are  admissible  are  not  completely  arbitrary,  but 
it  is  clear  from  the  experience  of  Karp,  Kotik,  and  Lurye  that  certain  singu¬ 
larities  may  be  permitted.  The  analogy  to  the  requirement  which  they  found 
necessary,  when  translated  from  a  dipole  distribution  to  a  source  distribution, 
is  the  requirement  that  the  potential  0’  remain  bounded  when  it  is  described 

0'  •  const  f  ff  G(x,y,s;  h,k,-f)  dV  ...(3.29) 

V  r 

where 

0  -  0'  +  cx  .  . . .(2.44) 

The  integration  V  is  over  the  volume  in  which  the  source  distribution  is 
non- zero  and  0'  must  be  bounded  only  outside  this  volume. 

Bxtctly  the  same  treatment  as  outlined  above  may  be  applied  to  distri¬ 
butions  of  dipole  moment  density  ji(h ,k ,-f )  rather  than  source  density.  This 
requires  only  that  we  substitute  equation  (2.33)  or  its  equivalent  in  terms 
of  an  integral  into  equation  (2.25)  to  get  the  resistance  in  terms  of  dipole 
density  rather  than  source  density.  Then  we  get  a  result  equivalent  to 
(3.28)  except  that  the  resistance  components  must  be  written  R(pr»j*#)  rather 
than  R(6r,6a).  In  such  a  case  the  restrictions  on  the  singularities  in 


44 


are  probably  given  by  equation  (3*15)  extended  to  include  the  case  k  i  0.  It 
la  also  posalble  to  use  assemblages  of  discrete  sources  or  discrete  dipoles  ^ 
in  place  of  the  densities  discussed  above. 

The  treatments  using  this  concept,  despite  their  simplicity,  have  been 
limited  mainly  to  work  which  was  started  by  G.  P*  Welnblum  at  the  David 
Taylor  Model  Basin,  and  carried  further  at  the  towing  tank  at  Stevens  Institute 
of  Technology  and  at  the  University  of  California.  Welnblum  described  his 
dipole  distribution  (which  he  did  not  initially  distinguish  from  the  offsets 
of  the  hull)  by  the  product  of  a  polynomial  in  the  variable  f  describing  the 

distance  below  the  surface  and  a  polynomial  in  the  variable  h  describing  the 

longitudinal  position.  After  converting  these  to  non-dimensional  form  he 
could  write  the  offset  7  the  form 

^-X(>)Z(j)  -  (1  -z:  «„£"  -Ebml")(l  -  etr)  ...(3.30) 

where  the  half-length  of  the  hull  is  1  so  that  -1  <  1  and  the  draft  of 

the  hull  is  1  so  that  0  <  ^  <  1.  The  signs  of  the  an  are  so  arranged  as 
to  make  them  describe  functions  symmetric  fore-and-aft,  while  the  bm  are  for 
anti-symmetric  functions.  The  quantity  e  <  1  is  a  positive  parameter  which 
is  used  to  permit  the  hull  to  have  a  flat  bottom  (if  e  <  1).  In  his  first 
report  |^5j  Welnblum  used  r  *  4  and  limited  the  calculations  to  n  *  2,3  ,4,6,8, 

10,12.  He  then  calculated  the  resistance  terms  which  enter  into  a  sum  of 

the  form  (3.25),  explicitly  differentiating  the  dimension  (or  dipole  moment  distri¬ 
bution)  to  convert  it  to  the  form  6  which  enters  into  that  equation.  He  also 
. 

used  another  distribution  in  the  direction  of  depth  which  gave  him  the 
equivalent  of  a  V-bottom. 
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In  Weinblum's  Initial  paper  the  integrals  corresponding  to  the  values 
of  n  and  r  mentioned  above  are  given.  Later  the  calculation  is  carried 

* 

through  to  the  point  where  optimum  forma  are  obtained  by  a  technique  similar 
to  that  described  earlier  [l6j.  These  forms  are  restricted  in  that  the  dipole 
moment  density  is  made  zero  at  the  ends  of  the  ship,  although  Karp,  Kotik 
and  Lurye  £6]  later  found  that  it  should  be  permitted  to  become  infinite  to 
minimize  resistance.  The  separability  of  the  functions  X(£  )  and  Z(^)  severely 
limits  the  forms  which  can  be  investigated.  Nevertheless,  even  these  forms 
show  the  possibility  of  significant  Improvement  over  conventional  fonts. 

For  example ,  Martin  and  White  £.12}  in  1961  selected  only  the  exponents 
n  *  2,3,4,  and  6,  added  the  constraint  that  the  value  of '7?  he  zero  at  bow 
and  stern,  constrained  the  Integral  ^  d  JF  »  constant,  and  then 

found  the  form  of  minimum  resistance  within  the  very  limited  family  of  shapes 
which  resulted.  Since  the  result  they  obtained  is  properly  Interpreted  as 
a  dipole  density  but  was  used  as  the  hull  offsets  in  building  a  model,  it 
is  clear  that  the  model  towed  was  probably  not  as  good  as  it  might  have  been 
had  the  streamlines  been  calculated  directly  from  the  dipole  density.  None¬ 
theless,  the  model,  which  was  optimized  for  a  Froude  number  of  0.5,  had  a 
wave -making  resistance  about  13  per  cent  less  than  that  of  the  model  of  a 
good  destroyer  hull  of  similar  dimensions  Intended  to  run  at  a  comparable 
speed.  The  success  of  this  limited  effort  shows  that  there  is  much  to  be 
gained  by  a  mathematical  approach  even  within  the  confined  of  conventional 
hull  shapes. 

In  an  effort  to  see  what  could  be  done  by  recognizing  that  the  form  of 
the  afterbody  of  a  ship  is  largely  determined  by  the  requirements  for  good 
flow  of  water  to  the  propellers,  J.  Richard  Gauthey  published  in  1961  the 
results  of  some  calculations  using  a  similar  method  specifying  the  afterbody 
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of  the  ship  completely  and  permitting  only  the  toward  half  to  vary  Jj.3^. 

Hia  results  bring  out  two  major  points;  (1)  the  form  which  resulted  was 
quite  different  from  the  forebody  which  had  been  designed  by  conventional 
methods  to  go  with  the  particular  afterbody;  and  (2)  small  changes  in  shape . 
can  have  large  effects,  Gauthey's  results,  while  limited  in  much  the  same 
fashion  as  those  of  Welnblum  and  of  Martin  and  White,  add  to  the  confidence 
that  much  can  be  done  to  find  Improved  hulls  of  conventional  form  by 
mathematical  methods, 

4.  CANCELING  WAVES  BY  ADDING  APPENDAGES  WHICH  PRODUCE 
SIMILAR  WAVES  OF  OPPOSITE  PHASE 

In  a  series  of  important  papers  starting  in  the  summer  of  1960,  Takao 
Inui  and  his  collaborators  at  the  University  of  Tokyo  described  the  results 
of  Intentionally  adding  large  bulbs  to  hulls  In  such  a  position  and  of  such 
a  size  as  to  produce  waves  about  of  the  same  form  but  of  opposite  phase  to 
those  produced  by  the  bow  and  stern  of  the  ship  The  ability  to 

do  this  might  be  inferred  from  equations  (3*5)  to  (3.8),  Those  equations 
show  that  a  hull  form  whose  equivalent  source  distribution  is  a  continuous 
function  of  position  along  the  hull  and  which  has  continuous  derivatives  of 
all  orders  produces  two  well-defined  wave  systems,  one  starting  at  the  bow 
and  the  other  at  the  stern.  Each  system  has  a  so-called  "sine  component" 
and  a  "cosine  component".  To  see  this,  we  may  write  the  first  term  of  (3.7). 

^ bow, sine  ■  4  S(0,©)  sln(K0x  sec2©  cos©)  cos sec2©  sin©>  d®  ...(3.31) 

c  -jt/2 

It  is  possible  to  assume  that  S(0,Q)  is  an  even  function  of  ©  for  a  ship 
which  is  syonetrical  on  the  centerline  plane,  as  nearly  all  ships  are.  Then 
since  sln(Koysec2©  sin©)  is  an  odd  function  of  ©,  it  follows  that 
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...(3.32) 


4  J*  8(0,0)  cos(K  x  sec29  cosO)  §in(K  y  aec^O  alnO)  dO  ■  0. 
c  -f/2 

We  aiy  now  add  equationa  (3.31)  and  (3.32)  and  gat 

X  bow, sine  -  4  $  S(0,0)  ain  [*K  (xcoaO  +  yalnO)  aac2©l  dO  ...(3.33) 

'  c  -rrfl  L°  J 

A  similar  manipulation  can  be  made  with  the  other  terma  to  provide  a  term  In 
coa  ^{^(xcoaOfyainO)  aec2oJ  and  a  pair  of  terma  with  x  — 9  x+L . 

Lunde  £  7]  has  shown,  that  quite  generally  the  vavea  produced  by  a  ahip  which 
ia  aymmetrical  on  ita  centerline  can  be  deacrlbed  at  pointa  a  long  way  from 
the  ahip  by  a  a urn  of  a  sine  component  and  a  coalne  component  as  follows: 

3 ahip  "  ~  ^2  S(h°’8)  •ln£Ko£<x-ho>co*«  +  yeinoj  sec20J.  dO  + 


c  -Mil 

rn 

+  A  S  C(h0, 

C  -*/2 


0)  cos  £k0[(x- 


hQ)cosO  +  yainO 


yalnoj 

jlnsJ  aec20^d 


••.(3.34) 


Where  the  functions  S(hQ>0)  and  C(ho>9)  are  the  results  of  combining  a  well- 
defined  bow  wave  and  a  well-defined  stern  wave,  they  are  not  likely  to  be  a 

form  which  is  easily  canceled  by  an  added  appendage  such  as  a  bulb.  On  the 

other  hand9  Inul  has  found  that  for  ships  of  certain  forms  moving  at  relatively 

low  speeds ,  the  wave  originating  at  the  bow  and  the  wave  originating  at  the 

stern  may  each  be  described  by  a  sine  component  whose  amplitude  function 
S(h,9)  is  positive  and  does  not  change  rapidly  with  increasing  9.  That  this 
is  likely  can  be  seen  from  equation  (3.6) t  which  clearly  has  as  its  dominant 
term  for  large  5^  (small  speed)  the  value 

S(h,0)  1  6(h)  ...(3.35) 

and  it  follow,  that  the  cosine  term  goes  to  aero: 

C(h,9)  — *  0. 

Kglarge 
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tkft  iUaeeift  ofc  IttUl1!  technique  is  based  bn  this  observation.  Ha  uses  the 
vell-kttoVti  tesUlt  that  the  waves  produced  by  a  doublet  can  be  described 
entirely  by  e  sine-wave  component.  The  amplitude  function  of  that  component, 
8(9),  is  not  really  very  close  in  form  to  that  of  a  ship's  bow  wave,  but  it 

is  opposite  in  sign.  By  matching  amplitudes  at  0  «  0,  and  counting  on  the 

factor  cos  0  in  equation  (3.8)  to  minimise  the  effects  of  the  differences 
for  large  0,  he  gets  virtually  complete  canceling  of  the  transverse  wave 
system  and  some  decrease  in  strength  of  the  diverging  system  of  waves.  Since 
most  of  the  energy  in  a  ship  wave  system  for  slow  and  intermediate  speed 

ships  is  in  the  transverse  waves,  this  results  in  a  striking  decrease  in 

wave-making  resistance  for  such  ships  whose  bow  waves  can  be  described  by  a 
sine  wave  system.  The  resistance  of  the  combination  of  the  sine  wave  bow 
waves  and  the  spherical  bulb  can  be  mitten  in  the  following  fashion  for  any 
given  speed; 


R  *  const 


1 T/2  (  x  . 

t  $  jAyOv©)  -  B  (hy.ey 


co.^9 


dO 


...(3.36) 


Here  ve  have  let  the  sine  wave  component  of  the  amplitude  function  for  the 
bow  be  written  A^(h^,Q),  where  h^  is  the  longitudinal  position  at  which  the 
bow  waves  originate.  The  function  -B(h^,9)  is  the  amplitude  function  of  the 
waves  produced  by  the  bulb,  and  h^  the  longitudinal  position  of  its  center, 
which  has  been  made  coincident  with  the  longitudinal  position  where  the  bow 
waves  originate.  The  amplitude  function  for  the  bulb  can  be  written 


4  2 

-B(hp,9)  -  -const  M  sec  9exp(-KQfsec  9) 


...(3.37) 


The  magnitude  of  the  constant  depends  on  the  convention  used  in  writing  source 
strength.  M  is  the  dipole  moment  of  the  spherical  bulb.  If  we  compare  (3.37) 
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with  (3.35)  where  we  have  for  a  very  deep  hull  «t  low  speed,  8(0)  ■  hy,0)  • 

const,  we  see  that  to  minimise  R  in  (3.36)  we  must  not  only  make  the  bulb 
such  a  else  that  the  quantities  and  B  have  roughly  the  same  absolute  value*  ' 
at  0  »  0,  but  we  must  make  the  product  sec^O  exp (-K^f sec^O)  as  nearly  constant 
as  possible  over  the  range  of  0  for  which  cos^O  has  a  significant  value.  The 
only  way  this  can  be  done  is  to  adjust  the  value  of  f,  the  depth  of  immersion 
of  the  center  of  the  bulb,  assuming  the  speed  and  so  to  be  fixed.  If  we 
set  K0f '*’'.55,  for  example,  the  value  of  B  will  be  about  the  same  for  0  -  35° 
as  for  8  -  0,  and  for  8  >  35°  the  magnitude  of  B  will  fall  off  rapidly.  In 

such  a  case  a  deep  ship  traveling  at  a  speed  of  32  ft/sec  (slightly  less 

2 

than  19  knots)  would  have  to  have  the  bulb  Inmersed  .55  x  32  /g  -  17.5  ft  to 
its  center. 

Although  Inul  has  only  used  spheres  to  provide  cancellation  of  the 
"sine  wave"  bow  waves  and  stern  waves  t  it  is  easy  to  show  that  a  source-sink 
pair,  the  source  forward  of  the  sink,  separated  by  a  little  less  than  a  half 
wave-length,  will  provide  a  wave  which  is  much  like  that  produced  by  a 
spherical  bulb.  To  see  this  consider  a  source-sink  pair,  the  source  forward 
of  the  sink  a  distance  6  and  at  the  same  depth  below  the  surface.  The  wave 
pattern  far  aft  of  this  source  which  is  produced  by  a  traveling  source  is 
described  by  the  following  equation: 

T/2  , 

\  ^  J  c(®)  cos  (K0P  sec  ...(3.38) 

^  -7T/2 

C(9)  ■  const  m  sec^O  exp(-KQf  sec^O)  ...(3.39) 

where 

p  *  x  cos  0  +  y  sin  9  ...(3.40) 

If  we  use  (3.38)  to  write  the  wave-pattern  produced  by  the  source-sink  pair, 
we  get 
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...(3.41) 


*S 


I rft  ,  rrfi  , 

|  C(8)  cot  (K0Pj  «ecz8)  <18  -  J  C(8)  cot(K0p2ttcz8)  d8 
-F/2  -T/2 


where 

p  ■  p  +  i  cot  0  and  p.  »  p  *  A  cot  8 
1  2  2 

We  nay  tlmpllfy  equation  (3.41)  by  uaing  (3.42)  and  ve  get 

. Wll 

^^-2  J  C(8)  tin(K0p  tecz8)  .inO^  tec  8)  d8 
-Iti  2  2 

We  may  now  combine  the  portions  of  the  amplitude  function  and  write 


...(3.42) 


...<3.43) 


J  ^(0)  ain(Kcp  sec20)  dQ  ...(3.44) 

1  -m 

where 

0^(0)  -  -2  const  m  sec^O  sin(KQ  k  8CC  ®)  exp(-KQf  *«c2®)  ..*(3.45) 

We  may  compare  the  amplitude  function  for  the  source-sink  pair  with  the 
amplitude  function  for  a  sphere 9  that  is ,  equation  (3*45)  with  (3.37)#  It 
is  clear  that  they  are  of  the  same  form  except  that  in  place  of  one  of  the 
terms  sec  0  for  the  sphere  there  is  the  term  sln(K^  sec  0)  for  the  source- 
sink  pair.  This  in  one  sense  is  an  advantage ,  since  it  gives  us  the  quantity 
4  to  manipulate  in  order  to  Improve  the  match  between  the  amplitude  function 
for  the  bulb  and  the  amplitude  function  for  the  ship's  bow  wave  or  stern 
wave.  On  the  other  hand,  it  makes  the  bulb's  performance  more  sensitive  to 
speed,  since  it  also  contains  the  factor  K^.  The  "bulb"  described  by  the 
source-sink  pair  in  this  instance  is  an  elongated  body  of  revolution,  rather 
like  a  blimp.  If  this  kind  of  bulb  can  be  used  for  a  propeller  shaft  housing 
or  faired  into  the  hull  it  may  be  an  improvement  over  a  sphere.  On  the  other 
hand,  if  it  merely  sticks  out  forward  like  a  ram,  it  is  no  improvement  at  all. 
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So  far  all  the  applications  of  Inui's  method  have  consisted  of  adding 
bulbs  to  cancel  the  "sine  wave"  component  of  the  bow  waves  and  stern  waves 
of  relatively  slow  ships.  Because  the  sphere 9  which  is  the  form  generated 
by  a  concentrated  dipole 9  produces  a  wave  form  which  is  of  form  similar  to 
these  "sine  wave"  components  of  the  bow  and  stern  waves 9  but  of  opposite 
phase  9  it  is  natural  to  exploit  its  use.  However  9  the  cosine  wave  component 
becomes  important  in  ships  of  higher  speed  than  cargo  liners  —  destroyers 9 
for  example.  Since  these  ships  combine  high  speed  with  shallow  draft 9  a 
spherical  bulb  is  an  inappropriate  appendage  to  add  in  any  case.  (At  high 
speed  it  would  have  to  be  deep  —  see  discussion  following  equation  (3 .3$.) 
Therefore  further  investigation  with  a  view  to  finding  some  other  way  of 
canceling  the  waves  of  cosine  form  is  in  order. 

A  first  observation  on  this  matter  is  that  the  waves  produced  by  an 
isolated  source  or  sink  are  waves  which  have  the  phase  of  a  cosine  component. 
The  waves  produced  by  a  sink  are  opposite  in  phase  to  those  of  a  sou  ce. 

Hence  if  a  ship  produces  waves  of  cosine  form9  the  proper  appendage  to  add9 
in  theory 9  is  an  isolated  sink.  Unfortunately  an  Isolated  sink  is  no  more 
a  real  entity  than  an  isolated  magnetic  pole;  it  must  always  have  associated 
with  it  a  source.  However 9  it  is  easy  to  see  that  if  the  hull  can  reasonably 
be  described  by  an  Isolated  source  forward  and  an  Isolated  sink  aft9  its 
wave-making  resistance  can  be  greatly  decreased  in  the  following  way:  add 
a  sink  of  strength  equal  to  half  the  forward  source  at  the  location  of  that 
source  and  add  another  source  of  strength  equal  to  the  source  at  a  distance 
about  a  half  wave-length  forward  of  it.  This  is  equivalent  to  splitting  the 
source  in  two  and  moving  half  of  it  a  half  wave-length  forward  of  its  initial 
position.  Then  the  resistance  of  the  forward  system  (disregarding  its  inter¬ 
action  with  the  stern  wave)  will  certainly  be  reduced  to  less  than  half  of 
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what  it  vat.  This  follows  from  equation  (2.17).  Suppose  ve  let  m  be  the 
strength  of  the  original  source  and  f  its  depth.  Then  its  resistance  will 


be 


p  -  2  2  -2KQf  f 

R  *  xorpK  m  a  °  j 

*  0  2  o 


...0.46) 


If  ve  cut  the  source  in  two,  end  place  half  of  it  a  distance  (h^-h^)  forward 
of  the  original  source  such  that  K  (h  -h  )  ■  “TT ,  then  the  resistance  of  the 
combination  will  be 

«  -  .2  I-2**'  7.-2K»ft<l+t>1'2  t*l«  £l  ♦  c[*Ut>1/2]ldt  (J 

The  quantity  in  brackets  for  0  t<  5/4  is  less  than  unity,  and  since  for 
any  reasonable  value  of  KQf  very  nearly  all  the  value  of  the  first  Integral 
will  be  attained  in  this  range,  R  will  be  much  less  than  half  its  original 
value.  The  same  result  follows  from  equation  (2.23)  without  calculation, 
since  the  interference  part  of  the  resistance  will  certainly  be  negative  in 
this  case. 

Unfortunately,  the  technique  outlined  here  is  not  practical  for  small, 
high  speed  ships  for  the  reason  that  a  half  wave-length  at  35  knots  is 
several  hundred  feet ,  and  it  is  not  reasonable  to  add  so  much  length  to  a 
small  ship.  It  therefore  appears  that  Inul's  method  or  any  obvious  variation 
is  not  likely  to  provide  the  answer  to  finding  very  high-speed  hull  forms  of 
low  resistance. 


47) 


5.  USE  OF  THE  GRADIENT  OF  THE  WAVE  RESISTANCE  FORMULA  TO  FIND 
LOCATIONS  WHERE  CHANGES  IN  HULL  FORM  GAN  BRING  IMPROVEMENT 
In  1936  Hogner  made  use  of  the  bilinear  character  of  the  wave  resistance 
formula  to  determine  where  changes  might  be  made  in  the  sectional  area  of  a 
ship  to  decrease  the  wave -making  resistance  £ll7.  In  his  one  published 
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article  on  the  subject,  he  derived  the  function  which  may  be  regarded  as  the 
gradient  of  the  functional  describing  the  resistance  and  promised  further 
development  of  the  idee  in  a  second  article  which  unfortunately  was  never  * 
published.  This  concept  leads  in  a  natural  way  to  the  method  of  steep 
descent  for  the  calculation  of  Improved  hull  forma.  It  is  therefore  possible 
that  failure  to  follow  up  on  Hogner's  paper  has  resulted  in  a  twenty -five 
year  delay  in  using  this  simple  method  to  find  fonts  of  decreased  wave-making 
resistance. 

To  understand  Hogner's  method,  we  may  consider  a  quadratic  functional 
with  a  symmetric  kernel  K(h,h'): 

R(p)  -  fj  p(h)K(h,h')^(h')  dh  dh' 

If  we  now  add  a  small  quantity  e  times  a  delta- function  to  yu(h)  at  the  point 
h^,  we  can  write 

R(ji  +  .SOij))  h)  +  *S(h)J  K(h,h')  (ji(h')  +  «6(h')j  dh  dh' 

-  Jj  fi{h)  K(h.h')  /i(h')  dh  dh'  + 

+  •  Jj»(h)K(h,h')  dh  + 

+  « jK(hpH')^t(h')  dh  + 

+  «2K(h1,h1)  .  ...(3.48) 

Because  of  the  symmetry  of  K(h,h'),  and  the  fact  that  e  is  a  small  number, 
we  can  write 

R(p  +  .$(1^))  -  R(ji)  «  2«  J  p(h)  KOi.hj)  dh  +  «2  ROtj.hj)  ...(3.49) 

But  the  difference  in  the  two  values  of  E  is  Just  the  change  which  results 
from  the  addition  of  a  small  increment  e6(h^)  to  the  given  function  ji(h) ;  and 
if  we  divide  by  c  we  have  what  may  reasonably  be  described  as  the  gradient 
of  the  functional  R  if  we  let  e— *0. 

}(hl)  -  grad  R  -  2  J/i (h)K(h .hj)  dh  ...(3.30) 
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This  may  be  related  to  the  problem  of  finding  an  improved  strut  of  infinite 
depth  by  identifying  the  kernel  K(h,h*)  with  the  kernel  Y0(Ko|h-h'  |  ),  which 
is  clearly  symmetrical ,  and  Identifying  the  function  ji(h)  with  the  dipole 
moment  distributed  on  the  centerline  plane.  It  follows  that  the  function^ 
gives  the  effect  of  a  change  In  dipole  moment  density  on  resistance  as  a 
function  of  the  coordinate  h,  the  longitudinal  position  on  the  hull.  To  make 
a  change  to  the  hull  form  which  will  decrease  the  resistance  while  holding 
the  total  dipole  moment  constant  (thereby  holding  the  volume  per  unit  depth 
nearly  constant)  ,  we  may  make  a  small  change  in  the  dipole  moment  density 
by  adding  a  quantity  f  J  and  then  subtracting  a  quantity  ((1/L)  J  J  dh: 


/>(0)(h)  +  ej  (h)  -  e 


dh 


...(3.51) 


It  is  evident  that 


//■ 

L 


(1) 


(h)  dh 


Sr 

L 


(0) 


(h)  dh  +  e 


jj|(h>  -  I  |  J(h')dhj  dh  -  J ^(0)<h>  dh 


...(3.52) 


Hence  we  have  found  the  condition  for  an  iteration  to  find  an  Improved  hull 

*  * 

form  while  holding  the  volume  constant.  It  is  clear  that  no  further  improve¬ 
ment  can  be  made  when  the  condition  is  reached  that  Jj(h)  ■  constant.  But 
this  is  just  the  condition  for  minimum  resistance  which  was  found  by  use  of 
the  calculus  of  variations  in  the  problem  as  solved  by  Karp,  Kotik ,  and  Lurye, 
namely  equation  (3.19)  in  different  notation: 


J^(h)Y0(K0|h-h'|)  dh  -X  ...(3.53) 


From  this  it  is  evident  that  the  end 
result  is  equivalent  in  the  two  approaches. 
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In  hia  paper  Hogner  proceeded  as  far  as  deriving  the  function  J  (h)  for 
several  Infinitely  deep  hull  forms  which  could  be  described  by  elementary  % 
functions,  but  failed  to  carry  the  calculation  further*  In  a  later  portion 
of  this  paper  the  calculation  of  an  Infinitely  deep  strut  of  reduced  wave- 
making  resistance  will  be  carried  through  by  this  method,  and  will  be  shown 
to  be  roughly  equivalent  to  the  result  of  using  the  calculus  of  variations 
on  the  same  problem*  An  extension  of  the  method  to  more  complicated  problems, 
with  hull  forms  varying  in  three  dimensions,  will  also  be  outlined* 
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IV.  A  Method  of  Steep  Descent  for  Developing  Improved  Hull  Forms 


It  appears  reasonable  to  conclude  from  the  work  done  to  date  In  finding  ^ 
hull  forms  of  decreased  wave-making  resistance  by  mathematical  methods  that 
significant  improvements  can  be  made  over  conventional  hull  forms.  However, 
all  the  applications  so  far  have  been  limited  in  their  scope.  The  classical 
calculus  of  variations  has  only  been  applied  to  the  case  of  infinitely  deep 
struts.  Moreover  it  seems  too  difficult  for  shapes  of  hull  which  vary  in 
three  dimensions.  The  Ritz  method  has  been  applied  only  in  finding  the  best 
of  a  very  limited  family  of  polynomials,  but  it  does  produce  improvement* 

The  method  of  Inui ,  which  consists  of  canceling  the  ship  waves  with  roughly 
similar  waves  of  opposite  phase  produced  by  a  submerged  sphere,  has  consider¬ 
able  potential  for  intermediate  speed  ships  which  can  handle  the  addition  of 
a  large  sphere,  but  it  does  not  seem  suitable  for  very  high  speed  ships  like 
destroyers.  Further,  it  has  the  appearance  of  "fixing*'  a  poor  hull  design 
rather  than  calculating  &  good  one.  It  is  clear  from  the  limitations  of  the 
work  already  done  that  there  is  room  for  a  method  which  can  work  improvements 
in  hull  forms  which  are  partly  constrained  in  their  shape  (for  example,  by 
the  shape  of  what  they  must  carry)  but  which  can  still  be  varied  in  some 
respects  so  as  to  decrease  their  wave-making  resistance. 

1.  REASONS  FOR  CHOICE  OF  A  METHOD  OF  STEEP  DESCENT, 

AND  A  DESCRIPTION  OF  THE  METHOD 

The  choice  of  a  method  should  not  be  determined  by  the  fact  that  there 
is  no  proof  that  the  form  found  by  it  will  have  an  absolute  minimum  of  resist¬ 
ance.  The  problem  appears  to  be  that  any  hull  shape  with  a  stationary  value 
of  wave-making  resistance  will  only  be  optimum  compared  with  other  shapes  of 
a  restricted  set.  Other  sets  may  contain  shapes  which  cause  even  less 
resistance  than  the  stationary  value  already  found.  Techniques  for  identifying 
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such  other  seta  are  therefore  moat  desirable.  It  follows  that  any  method 
vhidh  not  only  finds  Improved  forms  in  a  given  set  but  points  to  other  sets 
which  may  be  better  is  worth  investigating.  In  this  respect  the  so-called 
direct  methods  of  the  calculus  of  variations  seem  to  have  an  advantage  over 
the  classical  Indirect  methods.  They  are  attractive  for  hull  form  calcula¬ 
tions  for  several  other  reasons  as  well: 

(1)  The  description  of  a  ship's  hull  is  limited  in  practice  to  the  defini¬ 
tion  of  its  coordinates  at  a  finite  number  of  points.  Since  the  direct  methods 
of  the  calculus  of  variations  deal  in  general  with  a  finite  number  of  variables 
initially  (although  they  may  allow  the  number  ultimately  to  become  infinite), 
there  is  a  natural  relationship  between  the  normal  description  of  a  ship  and 
the  language  of  the  direct  methods. 

(2)  If  it  is  not  possible  to  permit  the  number  of  variables  to  become 

infinite  in  the  limit,  then  there  is  still  much  to  be  gained  by  dealing  with 

a  finite  number.  With  only  a  finite  number  the  question  of  convergence  is 

greatly  simplified.  If  the  result  obtained  with  a  finite  number  of  variables 
is  as  accurate  as  the  description  of  a  ship  can  be  in  practice,  then  the 
practical  limit  of  success  has  been  reached  anyway. 

(3)  The  mathematics  of  the  direct  methods  are  relatively  simple. 

(4)  The  direct  methods  are  susceptible  to  simple  applications  of 

constraints. 

Although  the  direct  methods  used  so  far  (by  Welnblum  and  others)  do  not 
have  all  these  advantages,  there  are  other  methods  available  by  no  means  as 
limited  as  those  already  applied.  Further,  the  mathematical  difficulties 
Involved  in  the  classical  calculus  of  variations,  when  applied  to  any  but  the 
simplest  ship  shapes,  make  some  alternative  essential  for  the  more  complicated 
problems  of  practical  ship  forms.  The  choice  therefore  should  be  made  among 
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the  available  direct  methods. 

The  direct  methods  are  not  in  fact  so  far  removed  in  concept  from  the 
classical  indirect  ones.  Richard  Courant  has  pointed  out  £l8^  that  there  is  ^ 
a  close  relationship  between  the  Culer  differential  expression  of  the  class* 
leal  calculus  of  variations  and  the  gradient  of  a  function  in  a  finite- 
dimensional  vector  space.  In  fact,  Euler's  differential  expression  may  be 
considered  the  gradient  of  a  functional  in  function  space.  In  view  of  this, 
a  reasonable  approach  to  the  ship  resistance  problem  is  to  write  the  resist¬ 
ance  as  a  function  of  a  finite  number  of  variables,  which  we  will  consider 
as  the  elements  of  a  vector  in  a  finite -dimensional  vector  space,  and  then 
examine  the  behavior  of  the  gradient  of  the  resistance  in  that  vector  space. 

If  we  start  with  the  variables  so  valued  as  to  represent  a  ship  and  vary  them 
along  a  trajectory  opposite  to  the  gradient  of  the  resistance  (or  as  closely 
to  that  direction  as  the  applied  constraints  permit),  then  we  should  produce 
an  improvement  in  the  hull  form  if  one  can  be  produced  by  continuous  varia¬ 
tion  of  its  defining  variables.  If,  in  addition,  we  follow  the  trajectory 
far  enough,  we  may  approach  a  stationary  value  of  the  resistance  functional. 

A  Theorem  on  the  Method  of  Steep  Descent 

This  discussion  leads  us  to  the  employment  of  a  method  of  steep  descent 

for  the  calculation  of  improved  hull  forms.  In  order  to  justify  this  for 

the  particular  application,  we  must  prove  the  following  theorem: 

Theorem:  If  C  is  a  non-negative  continuous  functional  of  the  elements  of  a 

finite-dimensional  vector  space  W,  and  the  elements  m  «  (m  )  of  the 

I  L  n 

space  are  differentiable  functions  of  a  parameter  t,  then  if  m  Is  unrestricted 
or  if  m  is  only  restricted  in  that  certain  elements  are  not  functions 

of  t  or  that  m  must  be  orthogonal  to  some  given  vector  a,  then  the  vector 
m(t)  will,  within  the  constraints  imposed,  have  some  direction  which  it  can 
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follow  which  will  provide  the  moat  rapid  decrease  or  least  rapid  Increase  In 

the  functional  C.  This  A  will  be  a  function  of  m.  If  the  direction  provides 

#  > 

a  decrease  and  the  vector  m  follows  the  trajectory  m<t)  ,  then  C  will  reach 
or  will  approach  asymptotically  a  stationary  value ,  and  this  stationary  value 
will  be  a  relative  minimum  value* 

Lemma  1:  Let  W  be  a  finite-dimensional  vector  space,  and  let  C(W)  be  a 

finite  non-negative  continuous  real  functional  of  the  elements  of  W.  Then 

if  we  write  grad  C  as  the  gradient  of  the  functional  C  with  respect  to  the 

elements  of  W,  and  if  all  the  elements  ra^  ,  m^,...,m^  of  the  vector  space  W 

are  differentiable  functions  of  a  parameter  t,  then  C(m  ,  m  ,  ...,m  )  goes 

12  n 

over  into  a  function  of  t,  and  we  can  write 
n 

C(t)  -JI  m.C  -  (dm,  grad  C) ,  where  dm  »  (m. ,  ra-,...,in  )  ...(4.1) 

1  “i  dt  dt  1  1  n 

Now  if  we  choose  dm  such  that  its  absolute  value  is  fixed  and  that  its  direc- 
dt 

tion  is  subject  to  certain  constraints,  but  within  these  constraints  is  such 
• 

as  to  give  C(t)  as  small  a  real  value  as  possible,  then  as  t  increases  C(t) 
will  ultimately  reach  or  asymptotically  approach  a  stationary  value, 
t 

C(t)  *  grad  C)dt 

0  dt 

Proof:  If  the  smallest  value  of  C(t)  which  we  can  obtain  subject  to  the  n 
constraints  of  the  problem  is  positive,  then  C(t)  is  a  (minimum)  stationary 

value  and  the  condition  of  the  lemma  is  satisfied.  If,  on  the  other  hand, 

•  • 

C(t)  is  negative,  then  C(t)  will  decrease  as  t  increases  until  C(t)  has  increased 

to  zero  or  until  C(t)  has  decreased  to  zero,  whichever  occurs  first;  and  in 

either  case,  since  both  C  and  C  are  continuous  functions  of  t,  a  stationary 

(minimum)  value  of  C  will  have  been  reached.  If  C  *  0,  this  is  a  stationary 

value  because  of  the  continuity  of  C  together  with  the  fact  that  C  is  non- 
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negative.  If,  instead,  C  asymptotically  approaches  zero,  then  C  will  approach 
a  stationary  value  as  closely  as  ve  please;  and  although  it  will  not  reach 
a  true  stationary  value,  it  will  be  such  for  all  engineering  purposes.  Notice 
that  this  proof  is  independent  of  whether  the  functional  C(W)  is  a  bilinear 
functional  'of  the  elements  of  W.  It  is  enough  that  it  be  non-negative .  Of 
course  the  functional  C(W)  is,  to  first  order,  a  quadratic  functional  of  a 
vector  p  £  W;  but  although  this  is  a  necessary  condition  if,  for  example,  a 

method  such  as  We  inblum^s  method  is  used,  it  is  not  necessary  here.  It  is 

\ 

evident  from  the  physical  conditions  of  the  problem  that  C(W)  cannot  become 
negative:  this  would  be  equivalent  to  a  negative  wave-making  resistance, 
and  in  turn  would  imply  the  addition  of  energy  to  our  ship  hull  from  a  pre¬ 
viously  undisturbed  ocean*  It  is  also  clear  from  the  equations  themselves. 
Hence  a  stationary  value  of  the  resistance  can  be  reached  by  allowing  C  to 
vary  as  a  function  of  increasing  t,  and  our  lemma  is  proved. 


Lenina  2:  If  certain  of  the  (m,  ,m  , . . .  ,m  ),  say  m,  ,...m  are  not  functions 

12  o  n  k 

of  t,  then  Lemma  1  holds  except  that  the  definitions  of  £  and  dm  are  altered 

dt 


to  eliminate  terms  with  these  subscripts. 

Proof:  This  follows  from  the  fact  that  »  0. 

The  Interpretation  of  this  Lemma  la  that  if,  for  example,  it  is  desired  to 
hold  constant,  to  first  order,  certain  dimensions  of  a  ship  (that  is,  to  hold 
constant  the  values  of  certain  sources  and  sinks)  then  this  can  be  done  and 
the  trajectory  will  provide  a  atationary  value  of  C(t)  subject  to  these  con¬ 
straints.  It  also  follows  that  if  we  write  m  in  terms  of  a  basis  in  which 


are  orthogonal  to  the  remainder,  then  we  may  find  C(t)  in  terms  of 

the  derivatives  of  the  remaining  basis  vectors  and  so  find  a  trajectory  which 

is  orthogonal  to  the  vector  components  m  , . . . ,m  • 

h  k 
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Lean*  3:  The  direction  of  m(t)  which  for  *  given  magnitude  of  m(t)  will  pro¬ 
duce  the  most  rapid  change  in  C(t)  is  that  m(t)  which  is  parallel  to  grad  C.  > 
Proof:  C(t)  *  (m(t) ,  grad  C).  If  we  set  m(t)  *  m(t)  lm(t)  I  ,  then 

|c(t)||  -//m(t)^  //grade//  (m(t) ,  grad  C)<  lw(t)|  ffgrA«l  C/l  ,  ...(4.2) 

FWi****  CK 

the  equal  sign  holding  only  if  m(t)  is  parallel  or  anti-parallel  to  grad  C. 

Lemma  4:  The  direction  of  m(t)  which  for  a  given  magnitude  of  m(t)  will  pro¬ 
duce  the  most  rapid  change  in  C(t)  if  m(t)  is  constrained  to  be  orthogonal  to 
some  vector 9  say  grad  V9  is  that  direction  which  lies  in  the  plane  (two- 
dimensional  vector  space)  which  is  spanned  by  grad  C  and  grad  V,  and  is 
orthogonal  to  grad  V* 

Proof:  By  the  Gram-Schmidt  orthogonalization  process  we  may  construct  an 

orthogonal  basis  for  the  two-dimensional  subspace  spanned  by  grad  V  and  grad 

C.  We  take  grad  V  as  the  first  of  our  basis  vectors.  Then  we  obtain  a 

second  basis  vector  orthogonal  to  grad  V  by  writing 

■  grad  C  -  (grad  C,  grad  V)  grad  V  •••(4.3) 

(grad  V,  grad  V) 

We  can  show  that  <L  is  orthogonal  to  grad  V  simply  by  writing  out  (<L»  grad  V). 

In  addition  we  can  show  that  <L  is  not  orthogonal  to  grad  C  by  writing 

(a,  grad  C)  -  (grad  C ,  grad  C)  -  (grad  Ct  grad  V )*  ...  (4 .4) 

(grad  V,  grad  V) 

If  grad  C  4  0,  which  we  showed  in  lemma  3  must  be  the  case  if  we  have  not 

already  found  a  stationary  value  of  C,  then 

(<i.  grad  C)  •  1  -  (grad  C,  grad  V)2 
(grad  C9  grad  C)  (grad  C 9  grad  C)  (grad  V 9  grad  V)  ...(4.5) 

But  the  second  terra  of  this  expression  is  always  non-negative  and  less  than 

unity  unless  grad  C  is  parallel  to  grad  V,  in  which  case  ctwlll  be  orthogonal 
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both  to  grad  V  (which  it  is  by  definition)  and  to  grad  C»  But  this  means 
that  unless  grad  C  and  grad  V  coincide  in  direction,  <vvill  not  be  orthogonal 
to  grad  C«  Now,  since  a*  is  both  in  the  two-dimensional  subspace  spanned  by 
grad  C  and  grad  V,  and  lu  also  orthogonal  to  grad  V,  It  must  be  parallel  to 
that  vector  m  which  is  orthogonal  to  grad  V  and  which  we  have  stated  will, 
of  all  vectors  of  that  magnitude  orthogonal  to  grad  V,  produce  the  largest 
value  of  jc|* 

That  this  is  the  case  can  be  shown  as  follows: 


Let 


grad  C  *  a (grad  V)  +  bet 

but  since  m  is  parallel  to  a,  we  can  write,  with  cl  -  cm, 
grad  C  -  a  (grad  V)  +  b(cra). 


But 

C  *  (grad  C,  m)  *  a  (grad  V,  m)  +  bc(m,m)  -  bc(m,m)  • 

Now  let  us  assume  that  there  is  some  other  vector  m*  of  the  same  magnitude 
as  m  so  that  (m9,m9)  -  (m,m)  but  such  that  |  (grad  C,  m9)j  >  |[grad  C,m)J  • 
Then 

■'  -  b9cL  +  p, 

where  (grad  V,p)  *  0  and  (p,<t)  *■  0 

•  #  • 
so  that  (p,m)  -  0.  It  follows  that  if  C9  is  the  value  of  G  corresponding  to 

this  vector  m9 ,  then 

C9  •  (grad  C,  m9)  -  (a  grad  V  +  bcii,  b9cm  +  p)  -  bb9c^(m,m). 

But  2  2 

(m9,m9)  *  b1  c  (m,m)  +  (p,  p)  -  (m,m)  by  hypothesis;  so 
b'c  -  J  1  -  (p,f)/(m  ,m)  <1  1.  Then  since 
C9  *  (b9c)bc(m,«n)  -  (b9c)  C  , 

it  follows  that 
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But  this  la  contrary  to  hypothesis,  ao  the  Lemma  la  proved. 


Proof  of  theorem:  By  Lemma  1  ve  showed  that  so  long  as  C(W)  is  a  finite, 

non -negative  continuous  real  functional  of  the  elements  of  W,  there  will  be 

a  stationary  relative  minimum  value  of  C(W) .  We  also  showed  that  this  value 

can  be  found  or  approached  asymptotically  by  allowing  each  of  the  elements 

of  W  to  become  a  function  of  some  parameter  t,  and  then  varying  the  elements 

as  a  function  of  t  so  that  C  decreases  as  rapidly  as  possible  for  a  given 

absolute  value  of  dm.  By  Lemma  2  we  showed  that  if  certain  elements  of  W 
dt 

were  held  constant,  the  theorem  still  held  subject  to  this  condition.  By 

Lemma  3  we  showed  that  the  direction  in  which  the  function  dm  should  be  varied 

dt 

to  produce  the  required  most  rapid  possible  decrease  in  C  subject  to  the 

restriction  of  Lemma  2  could  be  found.  By  Lemma  4  we  showed  that  even  if  the 

restriction  were  Imposed  that  the  trajectory  of  m(t)  be  such  that  some  other 

functional  of  the  elements  of  the  vector  space  be  held  constant,  a  direction 

of  dm  could  be  found  so  that  for  a  given  absolute  value  |  dm  j  the  most  rapid 
dt  I  dt  | 

possible  decrease  would  result  in  C(W) •  But  since  this  is  so,  our  Theorem 

is  proved. 

It  is  possible  to  derive  the  direction  which  dm  must  take  in  terms  of 

dt 

the  given  basis  vectors  if  two  constraints  of  the  sort  described  in  Lemma  4 

are  applied.  However,  it  results  in  complicated  algebra.  In  consequence  it 

is  simpler  to  rotate  coordinates  so  that  the  basis  for  the  vector  space  includes 

the  gradient  of  each  of  the  functionals  to  be  held  constant.  Then  the  rest 

of  the  basis  vectors  are  made  orthogonal  to  each  of  these  gradients,  and  the 

gradient  vectors  of  the  functionals  to  be  held  constant  are  disregarded  in 

the  calculation  of  the  gradient  of  C(W).  This  follows  because  in  the  exprea- 
•  # 

•ion  C  ■  (grad  C,  m)  there  is  no  contribution  from  any  of  the  elements  of 
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where 


grad  C  which  are  parallel  to  the  basis  vectors  grad  V^,  grad  V2>  •••» 
the  grad  are  the  gradients  of  the  functionals  which  are  to  be  held 
constant.  Observe  that  if  the  are  linear  functionals  of  the  then  their  ^ 
gradients  will  be  Independent  of  the  value  of  m  and  so  constant  in  direction 
throughout  the  vector  space  W. 

2.  PRINCIPLES  OF  APPLICATION  OF  THE  METHOD  OF  STEEP  DESCENT 
TO  THE  SHIP  WAVE  PROBLEM 

We  may  now  identify  the  variables  of  the  theorem  of  the  preceding  sec¬ 
tion  with  the  variables  of  the  problem  of  reducing  wave-making  resistance. 

The  variable  C  may  be  identified  either  with  the  drag  coefficient  Cw  or 
the  wave-making  resistance  R  of  a  ship,  since  either  is  a  non-negative  func¬ 
tional  of  the  elements  of  the  finite -dimensional  vector  space  in  which  the 
hull  is  defined.  We  will  choose  to  identify  C  with  R,  the  wave -making 
resistance.  The  elements  of  W  we  may  identify  with  the  sources  which  are 
used  to  generate  a  set  of  closed  streamlines  which  outline  the  hull.  The 
ship  Itself  then  is  described  by  a  vector  in  the  vector  3pace  W.  Since  any 
n  orthonormal  linear  combinations  of  the  sources  also  comprise  an  orthonormal 
basis  for  our  vector  space  W,  the  elements  of  such  a  basis  may  also  be 
identified  with  the  of  the  preceding  section.  It  will  be  shown  later  that 
linearly  Independent  combinations  of  the  sources (with  certain  restrictions) 
may  also  be  used  instead  of  orthonormal  combinations,  since  they  too  comprise 
a  basis  for  the  vector  space  W.  One  such  linearly  independent  basis  is  just 
the  summation 

J 

**7  ■  ^  ...(4,6) 

i*l 

for  the  case  of  a  set  of  sources  and  sinks  evenly  spaced  along  the  axis  of 
the  hull.  Notice  here  that  the  in  equation  (4.6)  are  just  the  sources 
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themselves,  end  the  *7j  are  to  be  identified  with  the  of  the  preceding 
section.  The  ^  may  be  interpreted  as  the  first-order  approximation  to  the 
offsets  of  the  hull*  The  vector  defining  the  ship  may  then  be  written  as 
f)  *  2*  Inu*  ^as  atlown»  this  is  a  Poor  approximation 

to  the  hull  form,  but  it  is  Just  the  approximation  which  is  obtained  by  a 
literal  interpretation  of  Michell's  classic  paper  £lj. 

Another  combination  is  a  linear  combination  of  linear  combinations!  and 
this  can  be  written 


n 

r 

k-1 


j 

Urn 

i-1 


ik 


. . .(4,7) 


Equation  (4.7)  can  be  interpreted  as  describing  (still  to  first  order)  the 
shape  of  the  ship's  hull  as  the  sum  of  a  set  of  n  other  first  order  approxi¬ 
mations  to  ship's  hulls.  Each  of  these  other  hulls  is  taken  as  a  basis  vector 
for  our  vector  space  W,  and  the  hull  described  by  the  ^  is  then  a  vector 
described  In  terms  of  these  other  ships  as  basis  vectors*  The  method  used 
by  Weinblum  can  also  be  interpreted  In  the  same  manner ,  except  that  in  his 
work  he  used  the  terms  of  a  polynomial  as  the  basis  vectors  for  the  descrip¬ 
tion  of  his  hull* 

The  fact  that  the  wave -making  resistance  R  is  non-negative  follows  from 
the  physical  considerations  of  the  problem*  For  the  case  of  an  ideal  fluid 
it  is  also  obvious  from  one  form  of  the  equation  for  the  wave-making  resistance 
of  a  continuous  distribution  of  x-directed  dipoles  in  a  half-strip* 


R  -  l6i^>K02  J  (I2  +  J2)  cosh2  u  du  ...(2.25) 

o 

Here  the  functions  I  and  J  are  real  so  the  integrand  and  the  Integral  are 
clearly  non-negative* 
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The  functionals  may  be  identified  with  any  of  the  quantities  which 
must  be  constrained  to  be  constant  in  order  to  make  the  statement  of  the 
wave -resistance  problem  meaningful.  For  example,  it  is  necessary  that  the 
sum  of  the  source  strengths  be  zero,  and  ve  may  identify  with  the  sum  of 
the  source  strengths.  A  possible  constraint  is  that  the  volume  of  the  hull 
be  constant ,  and  we  may  identify  the  functional  with  the  volume*  (In 
practice  it  will  probably  be  simpler  to  identify  it  instead  with  the  sum  of 
the  first  moments  of  the  source  strengths  in  the  direction  of  ship  motion, 
which  is  a  first  approximation  to  the  hull  volume.  This  sum  is  linear  in 
the  source  strengths  and  so  has  a  gradient  which  is  constant  everywhere.) 

In  addition  to  the  constraints  discussed  above,  there  is  one  other  set 
which  is  Inherent  in  the  way  the  problem  is  set  up.  This  is  the  location  of 
the  sources.  They  are  located  as  part  of  the  initial  information  in  the 
definition  of  the  problem,  and  the  number  of  linearly  independent  vectors  in 
our  finite-dimensional  vector  space  W  must  always  be  less  than  the  v  nnber  of 
source  locations.  This  follows  from  the  fact  that  the  sum  of  the  source 
strengths  must  be  zero  if  the  sources  are  to  generate  a  closed  hull  form. 
There  is,  however,  no  limitation  on  where  we  place  the  sources.  For  example, 
if  we  wish  to  find  the  effect  of  extending  a  hull  beyond  its  original  length 
we  may  postulate  a  source  of  zero  strength  at  some  point  forward  of  the  bow. 
Then  we  may  calculate  the  effect  on  the  ship  resistance  of  increasing  the 
strength  of  this  source  by  an  infinitesimal  amount.  If  this  results  in  a 
decrease  in  the  resistance,  we  may  conclude  that  a  bow  extending  out  to  this 
point  may  be  of  some  value.  It  is  this  ability  to  test  the  effect  of  sources 
placed  arbitrarily  which  permits  the  method  of  steep  descent  to  be  used  to 
explore  changes  to  the  hull  beyond  the  original  framework  of  its  definition. 
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Approach  to  the  Minimizing  Problem 

The  outline  of  the  method  to  be  followed  now  becomes  clear.  First  we  ^ 
choose  a  hull  form  as  the  starting  point  for  the  variation.  Then  we  describe 
it  in  terms  of  a  set  of  sources  and  sinks.  Next,  we  add  in  further  sources , 
but  of  aero  strength,  at  all  locations  where  we  think  it  might  be  reasonable 
to  have  them  in  the  final  ship  design.  In  our  next  step  we  select  those 
sources,  or  functionals  of  source  strength,  which  are  to  be  held  constant  in 
intensity.  Then  we  calculate  the  gradient  of  the  wave -making  resistance. 

Then  we  find  that  vector  direction  which,  subject  to  the  applied  constraints, 
is  most  nearly  parallel  to  the  negative  of  the  gradient.  Following  this 
direction  will  produce  the  most  rapid  decrease  in  the  value  of  the  wave-making 
resistance.  Then  we  vary  the  sources  and  sinks  in  strength  in  this  vector 
direction,  recalculating  the  direction  to  be  followed  at  such  short  intervals 
that  the  direction  changes  only  a  small  amount  between  calculations.  In 
this  manner  we  follow  a  trajectory  in  the  n-dimensional  vector  space  in  which 
the  ship's  hull  is  defined  until  we  reach  or  make  an  asymptotic  approach  to 
a  stationary  value  of  the  wave-making  resistance. 

This  simple  concept  need  not  be  confined  to  the  minimisation  of  wave- 
making  resistance.  There  is  no  reason  why  it  cannot  be  extended  to  minimiz¬ 
ing  the  sum  of  the  wave-making  resistance  together  with  the  other  components 
of  resistance,  subject  only  to  the  proviso  that  the  total  resistance  be 
written  as  a  non-negative  functional  of  the  source  strengths.  There  is  also 
no  reason  why  we  cannot  use  any  formulation  of  the  resistance,  whether  it 
provides  precisely  correct  total  resistance  or  not,  so  long  as  it  furnishes 
a  correct  or  nearly  ccrrect  description  of  the  direction  of  the  gradient  of 
the  resistance  in  the  n-dimensional  space  in  which  the  hull  shape  is  defined. 
This  means  that  it  may  be  feasible  to  use  a  partially  theoretical,  partially 
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empirical  formulation  of  the  resistance  —  and  even  change  it  as  the  vector 
defining  the  hull  form  follows  its  trajectory  through  its  n-dimenslonal 
vector  space  —  so  long  as  the  direction  of  the  gradient  is  correctly  described 


Demonstration  that  the  Resulting  Ship  Shape  Will  be  Smooth 
One  requirement  of  any  method  for  finding  hull  forms  of  decreased  wave¬ 
making  resistance  is  that  it  not  introduce  unwanted  discontinuities  into  the 
hull  surface.  This  is  necessary  to  prevent  increase  in  other  components  of 
the  resistance  than  the  wave-making  resistance.  The  description  of  a  ship 
in  terms  of  the  streamlines  of  the  potential  flow  arising  from  a  distribution 
of  sources  and  sinks  inside  the  closed  streamlines  insures  that  the  ship 
shape  will  be  fair  and  continuous.  This  follows  from  the  fact  that  the  poten¬ 
tial  arising  from  each  source  or  sink  is  a  solution  of  Laplace's  equation, 
and  the  potential  corresponding  to  the  entire  array  of  sources  and  sinks  which 
describes  the  flow  around  the  ship  is  simply  the  sum  of  the  potentials  of 
the  individual  sources  and  sinks,  hut  since  we  have  restricted  the  distribu¬ 
tion  of  sources  to  within  the  streamline  which  defines  the  ship's  hull 
surface  and  the  medium  outside  the  hull  has  a  constant  density,  the  solutions 
of  Laplace's  equation  outside  the  hull  are  everywhere  continuous.  From  this 
it  follows  that  the  sum  of  the  potentials  arising  from  a  finite  number  of 
sources  and  sinks  within  the  hull  boundary  is  also  continuous,  and  since 
Laplace '8  equation  is  linear,  the  sum  is  also  a  solution  of  Laplace's  equation. 
But  if  this  is  so,  then  each  of  its  derivatives  ie  also  a  solution  of  Laplace's 
equation  and  is  also  continuous,  since,  for  example, 


s72fMl  -ajy2*] 

[3xJ  J 


and  so  on  for  all  higher 
derivatives.  Since  the  direction  of  any  streamline  is  defined  by 


69 


where  i,  J,  and  k  are  unit  vectors  la 


the  directions  of  the  x,y,c  axes  respectively,  and  since  M,  and  9g  are 

3x  ^  9s 

continuous ,  it  is  clear  that  there  can  be  no  sudden  changes  in  direction  of 
streamlines  which  result  from  the  flow  from  a  finite  number  of  sources  and 
sinks  located  sway  from  the  streamlines. 

With  this  assurance  that  the  ship  forms  defined  by  a  finite  number  of 
sources  and  sinks  distributed  inside  the  streamline  which  bounds  the  hull 
will  be  continuous  and  fair,  we  may  proceed  with  confidence  to  use  the  dlstri 
butlon  of  a  finite  number  of  sources  and  sinks  to  define  our  ship  forms. 

They  will  at  least  be  reasonable  forms,  and  not  the  polygonal  approximations 
which  would  result  from  applying  the  same  procedure  to  the  offsets  of  the 
hull. 


Change  of  Basis 

It  was  pointed  out  earlier  that  if  the  quantities  to  be  held  constant 
are  linear  functionals  of  the  elements  of  W,  it  may  be  desirable  to  change 
the  basis  of  the  vector  space  W  so  that  it  includes  among  its  elements  the 
gradients  of  these  quantities.  Then  all  the  other  basis  elements  are  made 
orthogonal  to  these  gradients,  and  the  trajectory  through  W  of  the  vector 
describing  the  hull  is  calculated  disregarding  the  basis  elements  which  are 
parallel  to  gradients  of  the  quantities  to  be  held  constant.  That  is,  any 
linear  combination  whatever  of  the  new  basis  elements  other  then  those 
parallel  to  these  gradients  may  be  added  to  a  given  vector  without  changing 
the  quantities  to  be  held  constant.  This  will  continue  to  work  along  the 
entire  trajectory  if  the  directions  of  the  gradients  of  the  functionals  to 
be  held  constant  are  independent  of  position  in  the  vector  space.  That 


independence  is  assured  by  the  requirement  that  these  functionals  be  linear. 

The  process  of  changing  the  basis  may  sound  like  a  complicated  one*  In  ^ 

fact  it  is  not*  since  the  only  use  to  be  made  of  the  new  basis  vectors  is  to 

find  what  relation  must  exist  between  the  quantities  added  to  the  several 

sources  and  sinks  which  define  the  hull.  All  calculations  involving  the 

resistance  equations  themselves  are  made  in  terms  of  the  original  set  of 

sources  and  sinks,  changed  in  magnitude  but  not  in  position.  This  means  that 

the  calculations  of  the  resistance  integrals  need  be  done  only  once  y  and  may 

even  be  entered  into  the  problem  as  part  of  the  initial  information. 

Since  the  natural  description  of  the  hull  is  in  terms  of  the  intensities 

of  sources  at  predetermined  locations  9  the  change  in  basis  will  in  general 

make  the  new  basis  elements  linear  combinations  of  the  natural  ones.  We  may 

designate  the  old  basis  elements  as  a  , each  corresponding  to  unit 

1  n 

values  of  sources  with  strength  m1t...»m  .  Then  we  may  write  the  vector 

1  n 

which  defines  the  hull  in  terms  of  this  basis  as  m  *  m,o.+  ...  +  n  d.  We 

11  n  n 

may  write  the  new  basis  elements  as  Then  each  one  must  be  multi* 

plied  by  an  appropriate  scalar  80  ^at  vector  which  defines 

the  hull  will  become  ...  +rJ?nPn*  Then  it  follows  that  the  new  basis 

is  related  to  the  old  by  equations  of  the  form 

hm\civ\  ...(4.8) 

Since  in  the  original  basis  it  was  possible  to  write  the  gradient  of  the 
functional  C  as 

grad  C  *  dC  +  dC  <t-  +  ...  +  dC  <t- 
ol  dmj  *  <5nn  n 

we  may  define  a  gradient  in 

the  new  basis  as 
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grad  Cl  ■  bC  B,  +  be  B,  +  ...  +  BC  B  . 

I P  v  **?,  2  r% 

) 

If  the  basis  vectors  Pi>****Pn  *re  orthonormal  like  the  basis  vectors 
then  it  follows  from  the  results  of  ordinary  vector  calculus  that  the  direc¬ 
tion  of  the  gradient  vector  is  constant  regardless  of  what  orthonormal  basis 

is  used  to  write  it.  On  the  other  hand9  if  the  basis  vectors  B . ,B  are 

'  1  rn 

not  orthonormal  it  is  not  necessarily  the  case  that  the  gradient  vector 
grad  Cl  Is  in  the  some  direction  as  grad  Cl  .  In  particular,  if  the  basis 

vectors  B  f .  ..,B  are  not  orthogonal  but  only  linearly  independent,  then  it 
'1  '  n 

is  easy  to  find  examples  where  grad  cj^  is  not  parallel  to  grad  cj^.  never¬ 
theless,  for  a  given  magnitude  of  'J&y*  +  +  •  ^  a  vector 

motion  parallel  to  grad  cj^  will  provide  the  largest  change  in  C. 

To  find  a  relation  between  the  gradients  in  .the  two  bases  we  can  examine 
the  effect  of  making  a  small  change  d^  in  the  coefficient  of  ^  . 

h**?. \  dVt 

This  shows  that  a  small 

change  d*^>  in  the  coefficient  of  ^  corresponds  to  a  change  c^d1^?  in  the 

coefficient  of  each  of  the  cu  •  But  then  we  can  write 

k 

dc  ■  rg  <=lk"'71)  ■  u  |s  «lk>  *yv 

^  where  dC  is  the  small  change 

in  the  functional  C  resulting  from  the  small  change  d^  in  the  coefficient 
of  a*  .  But  by  definition, 
dC  -  dC  d 


and  by  equating  coefficients  of  dh^  in  this 
and  the  preceding  equation  we  have 
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...(4.9) 


2£ 

*7i 


W.  ukt  us*  of  the  hv  basis  in  the  following  fashion:  We  write 

th*  gradient  grad  Cl  by  use  of  equation  (4.9)  and  then  write  the  notion 

vector 


*7  "  y  lh+  >j?2p2  +  •”  +  ^n-jPn-j  +  °  ^n-j+l  +*“+  0  Pn 

where  the  lest  j  heals  vectors  ere  given  zero  coefficients  because  they  ere 
perallel  with  the  gradients  of  the  functionals  to  be  held  constant.  For  the 
non-zero  coefficients  we  substitute  the  corresponding  coefficients  of  grad  C 
and  then  we  ere  able  to  write 


(bc_  a t)p.  +  (ac_  At)a  +  ...  +  (&c_  a  t)Bn 
»%  2  d%-J  1 


...(4.10) 


The  At  is  a  negative  real  number  chosen  to  make  the  step  size  in  the  steep 
descent  correct.  It  must  be  negative ,  of  course ,  to  provide  a  decrf  ise  in 
C.  We  now  have  a  vector  A  in  a  subspace  VCW  which  is  orthogonal  to  the 
gradients  of  the  functionals  which  are  to  be  held  constant.  This  vector 
d^ia  parallel  to  that  portion  of  the  gradient  grad  cj^  which  is  Included 
in  V.  But  for  a  given  magnitude  of  this  is  the  direction  which  will 
produce  the  greatest  change  in  C9  since 


AC  -  (grad 

for  small  enough. 

Hence  we  have  found  the  direction  of  steepest  descent  within  the  given 
constraints. 

We  now  use  equation  (4.9)  to  apply  the  step  change  4  ^  to  the  source 

strengths  m, , ...»  m  .  We  write 
i  n 
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(dc_  At)^.  -  <*C_  At)£cik«,  -  £.  (M.  4  t  c  )«, 
0^  k  *  k  lk  * 


Then  we  sun  the  coefficients  for  ell  i  of  each  a,  end  edd  the  result  to  the 


coefficient  of  au  in  the  expression  m  ■  m  a*+  . . .+■  &>«..  +  m  *.  . 

k  11  n  n  n 

*".*-*2^  • 


•#•(4, 


This  Iteration  continues  with 

1-*  2,  2->3t...  until  a  minimum  is  reached  or  the  calculation  is  terminated 
for  some  other  reason* 


It  is  easy  to  show  that  the  direction  of  grad 


v 


different  from  that 


of  grad  Cl  for  a  particular  case.  Suppose  that  all  the  basis  vectors 

I*-' 

Pl,p2,ee«  are  identical  with  except  that  \  +  i 

0  <«  <,1.  Then  suppose  dC  «  bC  *  1,  all  other  dC  «  0*  Then  from  equa- 

dm  dm  dm 

k  k+1  J 

tlon  (4*9)  we  have  dC  ■  1  *  bC  »  1  +  a*  But  from  equation  (4.8)  the 

b%*i 

component  of  grad  C  in  the  direction  of  will  be  (l+e)e  ■  a  +  e2,  which 

will  be  much  less  than  1  for  c  small  enough*  But  for  any  a  ^  0  the  basis 
^|,***t^n  spans  the  vector  space  and  is  admissible 9  so  it  is  clear  that  we 
may  have  grad  Cl  not  parallel  to  grad  C 

Ip 

It  Is  clear  that  if  we  have  a  basis  ]b  such  as  the  one  just  described  and 

make  e  small  enough,  the  path  of  steep  descent  calculated  by  the  method 

described  above  will  be  essentially  orthogonal  to  the  direction  «P  until 

k+1 

a  stationary  value  of  C  with  respect  to  all  other  allowable  directions  in 
the  vector  space  is  approached.  Then  the  path  will  acquire  a  significant 
component  in  the  direction  and  proceed  more  directly  toward  a  true 

stationary  value  with  respect  to  all  motions.  It  is  therefore  evident  that 


ID 
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the  peth  of  steep  descent  is  not  independent  of  the  basis  vectors  which  we 
use  to  describe  the  hull.  Nonetheless,  it  does  not  appear  that  the  ultimate 
destination  of  the  path  is  necessarily  affected.  On  the  other  hand,  since  ' 
there  is  no  proof  available  that  the  solution  to  be  obtained  is  unique  except 
in  the  case  where  the  resistance  functional  is  bilinear,  it  may  well  be  that 
two  different  paths  will  result  in  two  separate  stationary  values,  just  as  a 
traveler  going  down  a  hill  may  end  tip  in  any  of  several  separate  valleys  in 
the  bottom. 

It  is  possible  to  select  a  set  of  linearly  Independent  basis  vectors 
for  the  description  of  the  hull  in  many  different  ways.  It  should  be  observed 
first  that  if  we  use  n  sources  to  describe  a  hull  there  are  really  only  n-1 
degrees  of  freedom  in  the  choice  of  their  values ,  since  their  sum  must  be 
zero.  If,  further,  we  wish  to  hold  constant  the  hull  volume  to  first  order, 
then  we  must  hold  constant  the  first  moment  of  the  source  strength  in  the 
x-directlon.  This  provides  another  constraint,  leaving  only  n-2  degrees  of 
freedom.  If  we  choose  n-2  basis  vectors  which  are  orthogonal  both  to  the 
gradient  of  the  total  source  strength  and  to  the  gradient  of  the  first  moment 
of  the  sources  in  the  x-directlon,  and  linearly  Independent  of  each  other,  we 
will  have  a  basis  for  a  subspace  V  included  in  W  in  which  every  vector  is 
orthogonal  to  theae  gradients.  It  will,  moreover,  contain  every  such  vector. 
It  follows  that  any  trajectory  through  this  subspace  will  leave  the  volume 
of  the  hull  unchanged,  at  least  to  first  order.  We  may  find  the  best  such 
trajectory  of  C(W)  through  the  subapace  V  by  taking  the  gradient  of  C  with 
respect  to  the  basis  elements  of  the  subspace  V  and  following  a  trajectory 
parallel  to  this  gradient. 

It  may  sound  as  if  it  is  a  difficult  problem  to  specify  a  basis  which 
is  orthogonal  to  the  gradient  of  the  sum  of  the  source  strengths  snd  to  the 
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gradient  of  the  first  moment  of  the  source  strengths  in  the  x-dlrection.  In 
practice  it  turns  out  to  be  quite  simple*  A  few  examples  will  make  this  clear* 
Suppose  we  have  a  hull  which  we  have  described  by  eight  sources,  equally 
spaced  along  its  centerline  plane.  We  may  describe  their  strength  by  multiples 
of  eight  unit  sources  (o^,»*#,a)  which  comprise  a  basis  for  the  description 
of  any  ship  which  can  be  described  in  terms  of  sources  placed  at  these  loca¬ 
tions*  We  may  have,  for  example, 

2q*  4au  2a#  CL/  -a,  -2<v  -4cu  -2a^  • 
12345  6  7  8 

Since  the  sum  of  the  source  strengths  is  aero,  this  set  generates  a  closed 

set  of  streamlines  which  could  be  a  ship*  The  antl-symnetrlcal  form  makes 

this  generate  a  set  of  streamlines  which  are  symmetrical  fore-and-aft* 

We  may  provide  a  basis  (B  ,  ...,B  ),  such  that  any  vector  described  in 

'1  '6 

terms  of  its  elements  can  be  added  to  our  given  hull  without  changing  either 
the  sum  of  its  elements  (which  is  aero)  or  its  first  moment,  which  is  roughly 
proportional  to  the  volume  of  the  hull*  The  basis  might  be  the  following: 

^  *0^/2  -  aj  4-  o^/2  4*  Oo.^  4-  Oqj  4*  Oo£  +  +  Oog 

^2  —  4-  a^/2  -  o^4*o^/2  4“  Ool^  4*  Oo^  4*  Oo^  +  Oog 

Ba  ■  0(k  4*  0(t  4-  o-  /2  -  S;  +  Ol-/2  +  Oqu  4"  Odu  4*  Ool* 

'  ^  l  2  3  45  678 

B  m  Oo^  4-  0<l»  +  Oo-'  4”  O-* /2  -  au  +  <x//2  +  0<L/  +  Oou 
r4  1  2  3  4  5  6  7  8 

B.  ■  Ov  +  Oou  4*  Oa-r  4*  Oa#  4*  a~/2  •-  /2  +  Ocl* 

3  1  2  3  4  5  67  8 

B  *  Oou  4*  Oau  4*  0<x#  4*  0<L/  4“  Oo**  4*  cl»/2  -  ou  +  ou/2 

r6  1  2  3  4  5  6  7  8 

The  six  vectors  B_,...,  B  are  linearly  independent  of  each  other.  When  any 
rl  r6 

multiple  of  any  of  them  is  added  to  the  set  of  vectors  which  describe  the 
ship,  neither  the  sum  of  the  source  strengths  nor  their  first  moment  is  changed. 
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Hence  these  six  vectors  meet  the  requirements  we  have  set* 

We  may  also  define  vectors  which  hold  constant  othor  features:  for 
example,  we  may  require  that  an  appendage  which  can  only  be  produced  by  a 
set  of  three  sources  in  a  particular  magnitude  relation  to  one  another  retain 
its  shape  although  not  its  site.  To  do  this  we  may  make  one  of  our  basis 
vectors,  or  one  term  of  a  basis  vector,  the  sum  of  the  three  sources  each 
multiplied  by  its  appropriate  magnitude* 

We  can  generalize  our  basis  vectors  which  maintain  the  volume  of  a  ship 
constant  (to  first  order)  to  the  case  where  the  ship  is  described  by  a  number 
of  rows  of  sources  distributed  over  the  centerline  plane  of  the  ship.  A 
single  basis  vector  might  be  described  by  the  following  array  of  numbers 
(the  multiplied  basis  vectors  in  terms  of  which  the  hull  is  described  are 
omitted) * 

1  -1  0000000... 

P  “ 

1  0  00  -1  10000... 

0  00000000... 


We  provide  n-2  vectors  of  this  sort.  To  do  this  easily  we  may  hold  the  pair 
of  numbers  in  the  upper  left-hand  corner  constant  and  associate  with  it  in 
turn  n-2  pairs  like  the  one  shown  in  the  second  row,  each  with  its  axis  parallel 
to  the  direction  of  motion  of  the  ship  and  with  the  order  of  signs  reversed 
from  the  order  of  the  reference  pair  in  the  upper  left-hand  corner.  Observe 
that  some  of  the  additional  pairs  will  have  one  member  in  each  of  two  rows. 

Also  ve  must  include  one  vector  which  looks  like  this  (we  will  call  it  JJ  ): 

-  1  -1  1  0  0  0  0  .  .  .  1  -2  1  0  0  0  ... 

0  000000...  000000... 

0  000000...  000000... 
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If  ve  wish  to  hold  the  sectional  eree  curve  constant  ve  nay  use  a  set  of 
basis  vectors  like  these: 


1  -1  0  0  0  0 

p  ■  -1  1  0  0  0  0 

0  0  0  0  0  0 


1  -1  0  0  0  0 

p  ■  0  0  0  0  0-  0 

-1  1.0  0  0  0 


^-01-1000 


-  0  1-1  0  0 


0 


0-11000 


0  0  1  0  0  0 


0  0  0  0  0  0 


0  -1  0  0  0  0 


Observe  that  if  there  are  m  columns  in  the  last  array  that  there  will 
be  a  total  of  n-m-1  linearly  independent  basis  vectors*  Hence  if  the  sec¬ 
tional  area  curve  is  kept  constant  there  will  be  fewer  degrees  of  freedom 
than  if  only  the  hull  volume  is  held  constant.  It  is ,  of  course ,  possible 
to  provide  sets  of  vectors  which  hold  the  volume  constant  on  each  waterline 
or  which  observe  other  constraints • 


Iteration  Sequence 

Once  we  have  described  the  ship's  hull  in  terms  of  sources  and  sinks  ve 
may  cow  calculate  the  resistance  by  one  of  the  formulae  developed  in  section 
II.  Then  we  determine  what  is  to  be  held  constant  and  develop  a  set  of  basis 
vectors  for  the  description  of  changes  in  the  ship's  hull  with  this  in  mind, 
as  described  in  the  preceding  paragraphs.  Then,  assuming  We  have  decided  to 
identify  the  wave-making  resistance  with  the  variable  C,  we  calculate  the 
partial  derivatives bR.  Then  we  calculate  the  partial  derivatives  of  R  with 

&BI 

respect  to  the  magnitude  ^ ? '  of  the  new  vectors  (except  those  identified 

with  gradients  of  vectors  to  be  held  constant)  using  equation  (4.9), 
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substituting  R  for  C.  If  there  «re  2  constraints  wo  will  have  n-2  such 
partial  derivatives.  Finally  we  can  write 

-At  21  pt  ■  4t  Z.  t-C  .o^  ,  &t<0.  ...(4.12) 

i  i  k 

v^(J+l)  .  rr?U)  +4^J)  ...(4.13) 

This  expresses  in  vector  form  the  same  result  as  is  given  in  equation  (4.11). 

The  simple  iterative  formula  (4.13)  serves  to  lead  us  to  the  desired  form  of 
relative  minimum  ^resistance  by  a  path  of  steep  descent  down  f rom  y  ^ .  It 
is  necessary  to  insure  in  the  Iteration  that  the  steps  are  kept  short  enough 
so  that  the  minimum  point  is  not  overshot;  or  that  the  change  in  direction 
of  the  actual  gradient  over  the  length  of  the  step  is  not  so  great  that  a 
step  which  makes  for  a  decrease  in  resistance  In  its  initial  portion  results 
finally  in  an  increase.  The  simplest  solution  to  such  difficulties  is  to 
proceed  until  such  are  found ,  and  then  return  to  the  preceding  step  and  try 
again  with  the  step  length  At  cut  in  half.  The  magnitude  of  succes  >ive 
values  of  the  double  sum  in  (4.12)  can  be  used  as  an  indication  of  whether 
progress  is  being  made  toward  a  stationary  value  of  R.  If,  for  exanq>le,  the 
direction  of  the  vector  changes  but  its  magnitude  does  not  ultimately 

decrease,  then  it  is  conceivable  that  the  iteration  is  resulting  in  a  circular 

m * 

f'  ' 

path  instead  of  following  the  path  of  steepest  descent;  but  if  this  occurs, 
then  the  size  of  the  steps  may  be  cut  in  half  in  order  to  make  the  actual 
path  followed  stay  closer  to  the  path  of  steepest  descent. 

For  the  case  where  all  derivatives  higher  than  the  second  are  zero  it 
is  possible  to  calculate  a  value  of  ^t  which  will  provide  the  largest  possible 
decrease  in  the  functional  C.  Since  the  resistance  functional  is  bilinear 
this  will  be  the  case  in  the  wave-making  resistance  problem.  Suppose  that 
we  can  write  AC  in  the  form  of  a  Taylor  Series  of  only  two  terms: 
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ac-  z  ac  (&t  be  \  +  i  r  x  ~fr2c  'bc_\ut  be  \. 

r-l®<vpry  T>yxJ  2l  .-1  r-1  byt  ) 

Since  At  <0,  we  nay  differentiate  with  respect  to  At  and  set  the  result 
equal  to  zero  as  follows: 

j(ac)  -  z.(bc_x  •  2  |  At  I  x  tL  B2c  be  be  -  o 

<1(40  r»l \®yry  2! I  I  *-l  r-1  d<?£ 

and  so  ve  have  for  the  optimum  size  of  the  step  4  t 

2 


optimum 


n  /  x  , 

2L  /JC  V 
r-i 


n  n 
21  21 

s-1  r 


dZC  dC  *hc 


1 


A  short  discussion  of  the  computational  techniques  Involved  can  be  found 
in  a  chapter  by  Charles  B.  Tompkins  in  ,fModern  Mathematics  for  the  Engineer" 

M- 

It  is  also  possible  to  use  Instead  of  the  rotation  of  coordinates  the 
constraint  of  equation  (4,3)  to  provide  a  trajectory  which  is  orthogonal  to 
the  gradient  of  a  function  which  is  to  be  held  constant.  This  has  the  dis¬ 
advantage  that  it  is  really  simple  only  if  there  is  a  single  constraint  -- 
not  a  common  situation. 

It  may  be  asked,  "Why  not  just  differentiate  the  resistance  with  respect 
to  each  of  the  source  strengths  in  turn  and  set  each  to  zero,  then  solve  the 
set  of  resulting  linear  equations?"  The  answer  is  that  it  may  be  undesirable 
to  permit  the  hull  shape  to  change  all  the  way  to  the  point  of  minimum 
resistance.  Since  the  gradient  will  decrease  as  the  relative  minimum  of 
resistance  is  approached,  a  larger  and  larger  change  in  hull  shape  will  be 
required  to  obtain  a  given  decrease  in  resistance  as  the  iteration  proceeds. 
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Presumably  the  iteration  will  start  with  a  conventional  hull  form  and  show 
how  it  may  be  changed  to  decrease  its  resistance.  There  will  probably  be 
some  hull  form  found  whose  wave-making  resistance  is  less  than  that  of  the 
original  and  acceptable  from  other  standpoints,  and  it  may  well  be  that  iteration 
beyond  that  point  will  produce  a  form  which  is  unacceptable.  It  is  also 
true  that  the  path  of  steep  descent  used  is  not  required  to  be  the  path  of 
steepest  descent.  That  is,  certain  of  the  basis  vectors  may  be  varied  less 
than  required  by  equation  (4.12),  or  not  varied  at  all,  so  long  as  the 
result  of  the  variation  is  a  decrease  in  total  resistance.  This  may  well 
produce  an  acceptable  form  in  cases  where  the  strict  method  of  steepest 
descent  produces  an  unacceptable  one. 

3.  THE  RESULTS  OF  A  SIMPLE  APPLICATION  OF  THE  METHOD  OF  STEEP  DESCENT  TO  THE 
PROBLEM  OF  IMPROVING  AN  INFINITELY  DEEP  PRISMATIC  HULL 

In  order  to  try  out  the  method  of  steep  descent,  a  sample  calculation 
was  made  using  a  slide  rule.  This  calculation  started  with  a  hull  which 
could  be  described  by  a  uniform  distribution  of  dipoles  with  axes  oriented 
in  the  direction  of  motion  over  a  centerline  plane  of  unit  length  and  infinite 
depth.  This  is  equivalent  to  a  two-dimensional  source  near  the  forward  end 
of  a  hull  and  a  two-dimensional  sink  near  the  after  end.  A  Froude  number  of 
0.316  was  assumed  for  the  calculation.  The  calculation  was  constrained  to 
hold  the  total  dipole  moment  constant.  After  six  iterations  the  shape  of 
the  dipole  distribution  approached  that  which  Karp,  Kotik,  and  Lurye  ^6} 
calculated  for  the  same  problem  with  Froude  number  0.38  using  the  classical 
calculus  of  variations,  and  differed  from  it  in  the  direction  to  be  expected 
from  extrapolation  from  their  results  for  higher  Froude  numbers. 

This  problem  has  only  one  constraint,  so  the  method  of  equation  (4.3) 
can  be  used  to  handle  it.  That  is,  with  the  resistance  R  and  the  sum  of  the 
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dipole  aonent  V,  a  vector  a.  parallel  to  the  dealred  direction  of  change  a(t) 

could  be  calculated  by  , 

a >•  grad  R  -  (grad  R.  grad  V)  grad  V  ...(4.3a) 

(grad  V,  grad  V) 

The  equation  uaed  for  the  reeietance  vaa  equation  (3.17), 

.,1/2  1/2  .  |  v  V 

/>*•  ./„  ,/2  p<j > /■<  p  V'Jt  1 • 

A  foci  symmetrical  fore-and-aft  was  assumed  at  the  start,  and  since  this 
provided  symmetrical  changes  the  symmetry  persisted  through  the  calculation. 

This  nearly  halved  the  amount  of  calculation  required  over  what  an  asyissetrlcal 
distribution  would  have  required.  The  length  of  the  dipole  distribution  was 
divided  into  ten  intervals,  and  with  one  exception  the  dipole  moment  over 
the  whole  interval  and  the  value  of  Y^  was  assumed  to  be  that  for  the  argument 
at  the  mid-point  of  the  interval.  The  exception  was  the  case  where  this  pro¬ 
duced  the  result  Yo(0)  ■  -®*.  Here  it  was  necessary  to  use  a  more  precise 

calculation,  since  the  function  YQ(x)  has  a  singularity  at  x  »  0  but  the 
x 

integral  J  YQ(x)dx  remains  bounded  for  all  x>0.  The  expression  for  the 
o 

Integral  is  available  in  terms  of  other  tabulated  functions  as  follows: 
x 


S  T0(n)d*  "  *  YQ(x)  +  2  x  j" Yj(x)  H0(x)  -  Y0(x) 


HjCx)! 


...(4.14) 


Here  Y^(x)  is  a  Bessel's  function  of  the  second  kind,  and  Hq(x)  and  H^(x)  are 
Henkel  functions. 

The  length  of  each  Interval  of  the  distribution  was  taken  as  A'g  and  the 
Intervals  were  numbered  from  1  to  10  with  midpoints  \  ,J  ■  l,2t..*,10.  The 
dipole  density  at  midpoint  was  written  M  •  Then  the  approximation  to  the 


J 


resistance  was  written  as 
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*  -  -8» 2  p  K02£ 


\  [5i  “‘Ki  ’•‘"•I  5  v + hj1  £  J*ot*W><*j4ii 


.<4  US) 


Then  by  differentiation,  and  changing  tha  limits  on  tha  Integral , 

g  ■  -«*’  f  *.2U  »!  t0»«,|5  *  “j  f  J\«4  *|>  *145, 

j  L  ...(4.16) 


It  vaa  then  faaaible  to  write 
grad  R 


(&»&»•..»%  l 

»*1  «*2  "nj 


...(4.16) 

...(4.17) 


81nce  (tin  aquation  (4.3a)  la  linear  In  grad  R,  and  we  are  Intereated  only 
In  Its  direction,  we  nay  disregard  the  magnitude  of  the  constants  In  equation 
(4.16)  in  calculating  a.  Further,  since  grad  V  enters  both  denominator  and 
numerator  to  the  asms  power  In  (4.3a)  we  may  multiply  grad  V  by  any  convenient 
multiplier.  If  we  take -  constant  for  all  j,  then  we  may  write  (grad  V) ^  ■ 
■  1  for  all  j.  Here  (grad  V) ^  la  the  J**1  component  of  grad  V.  With  this 
set  of  simplifications  we  can  now  write  the  equation  for  o* 


Observe  that 


■  2 

(grad  V,  grad  V)  -'  TL  (grad  V),  -  n,  and 

j-l  J 


n  n 

R ,  grad  V)  ■  £  bR  .■£  bR  . 

j-l  2Mj*  1  j-l  SMj 


(grad 
It  follows  that 


a»j  -  bR  -  1  ZL  bR 


where 


dHj  n  j-l  dMj 

ll  1  "J 

and  by  disregarding  the  multiplicative  constants  In  (4.16)  we  get 


...(4.18) 


...(4.19) 
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...(4.20) 


i  X-  M' 

J  5  n  J-l 

J  J 


dropping  the  constants. 


whore  1  and  R9  are  the  results  obtained  by 


The  procedure  used  is  to  calculate  the  bR9  from  equation  (4.16) ,  then 


the  o^1  fron  (4.20)  ,  and  then  find  a  new  value  for  the  m ^  by  taking 


...(4.21) 


Since  the  multiplicative  constants  were  negative  in  (4.16),  we  require 
At>0  for  this  iteration.  It  continues,  of  course,  with  0-*k,  l-*k+l  in  (4.21). 
The  limitation  on  selecting  the  value  of  A  t  was  to  make  sure  that  it  was  not 
so  big  as  to  make  any  of  the  m^  become  negative.  Otherwise,  the  slse  was  an 
estimate  by  the  calculator  that  it  might  produce  rapid  convergence  of  the 
calculation.  As  a  measure  of  the  progress,  it  is  possible  to  use  the  magnitude 
of  jaJgj-  >/(aJ,cJ  ) .  In  the  test  calculation  the  values  of  |  a*  |  went  in  the 
ratio  .82  :  3.74  :  2.02:  1.04:  .69:  .41:  .24.  The  initial  Increase  was  the 
result  of  the  first  step  being  too  large*  Later  decreases  resulted  from  using 
smaller  steps  —  from  a  little  experience  on  the  part  of  the  calculator.  The 
last  value  in  the  sequence  is  the  value  of  |  cuv  |  computed  from  the  values  of 
(These  are  the  final  values  of  the  calculation.)  This  means  that  the 
dipole  density  distribution  reached  in  the  particular  calculation  was  not  an 
actual  stationary  value,  but  was  certainly  a  much  better  shape  than  the  one 
with  which  the  calculation  was  started.  In  order  to  see  how  it  compares  with 
others  developed  by  the  calculus  of  variations,  it  is  plotted  together  with 
two  distributions  for  higher  values  of  the  Froude  number  calculated  by  Karp, 
Kotik,  and  Lurye  |j>J. 
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Figure  IV-1.  Dipole  Density  for  a  Strut. 

4.  USE  OF  THE  SOURCE  DISTRIBUTION  CALCULATED  BY  THE  METHOD 
OF  HESS  AND  SMITH  TO  FIND  IMPROVED  HULL  FORMS 
It  was  pointed  out  in  part  II  of  this  paper  that  one  of  the  more  difficult 
parts  of  the  calculation  of  hull  forms  of  reduced  wave-making  resistance  is 
to  find  a  source  distribution  equal  to  a  given  hull.  The  most  useful  such 
distribution  is  one  in  which  the  sources  are  placed  on  the  centerline  plane 
of  the  hull.  However ,  there  is  no  simple  general  method  for  calculating  a 
centerline  plane  distribution  corresponding  to  a  given  hull,  although  there 
are  many  special  cases  available.  On  the  other  hand,  there  is  a  general 
method  available  for  describing  an  arbitrary  hull  form  by  sources  placed  on 
its  surface,  and  this  method  will  provide  a  basis  (both  figuratively  and 
literally)  for  exploring  in  a  general  fashion  the  desirability  of  appendages 
to  add  to  a  given  hull.  Such  appendages  are  not  limited  to  spheres,  but  can 
be  virtually  any  shape  which  can  be  described  by  sources  and  sinks. 

The  method  of  Hess  and  Smith  £9],  which  Is  described  in  part  II  of  this 
paper,  provides  a  description  of  an  arbitrary  hull  in  terms  of  the  strengths 
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of  sources  and  sinks  distributed  in  a  regular  network  over  the  surface  of  a 
hull.  It  is  possible  to  describe  the  theoretical  wave-making  resistance  of 
such  a  set  of  sources  and  sinks  by  equation  (2.17).  (If  a  model  of  the  ship 
has  been  towed,  it  is  probable  that  the  theoretical  resistance  curves  can 
be  brought  into  precise  conformance  with  the  experimental  curves  by  use  of 
Inul's  semi -empirical  parameters.)  With  such  a  description  it  is  possible 
to  find  the  effect  of  additions  to  the  hull ,  whether  Incorporated  into  the 
hull  shape  or  Installed  as  exterior  appendages. 

To  explore  the  effect  of  appendages ,  all  that  need  be  done  is  to  add  to 
the  set  of  sources  describing  the  hull  a  set  of  sources  of  zero  strength 
distributed  in  the  locations  in  which  it  is  desired  to  explore  the  effect  of 
additions.  A  spherical  appendage  can  be  described  by  a  source-sink  pair, 
the  source  a  very  short  distance  forward  of  the  sink.  Other  appendages  may 
be  described  either  by  sets  of  simple ,  Independent  sources  placed  in  the  proper 
locations,  or  sets  of  sources  linked  by  definite  amplitude  relationship, 
depending  on  the  restrictions  to  be  placed  on  the  form  of  the  appendage.  In 
either  case,  the  method  of  steep  descent  is  then  applied,  permitting  only  the 
sources  being  investigated  to  vary.  If  the  sources  are  linked  by  a  definite 
amplitude  relationship,  then  it  will  be  necessary  to  apply  equations  (4,9) 
and  (4.10)  to  find  the  particular  basis  vector  for  the  change;  otherwise  they 
are  handled  without  rotation  of  coordinates ,  which  corresponds  to  c  ■  1 
for  k  •  i,  -  0  for  k  4  i  in  equation  (4.9).  In  any  case,  we  proceed  to 
find  the  values  of  dR  using  (2.17),  (4.9)  and  (4.10)  and  then  use  the 


method  of  steepest  descent  in  which  we  vary  only  the  sources  and  sinks  which 
are  added.  It  will,  of  course,  be  necessary  to  constrain  these  added  sources 
and  sinks  to  have  a  total  source  strength  of  zero  in  order  to  retain  closure 
of  the  hull  form.  It  will  also  be  necessary  to  Insure  that  the  variation  is 
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such  that  at  no  point  is  a  net  negative  change  made  in  the  firat  moment  of 
the  aourcea  in  the  direction  of  motion.  Thia  would  probably  result  in  the 
envelope  of  ihe  main  hull  moving  inside  the  set  of  sources  originally  used  > 
to  describe  it,  and  leave  them  to  describe  an  isolated  appendage  Instead. 

Also,  an  Isolated  source-sink  pair  with  the  sink  ahead  of  the  source  is  a 
meaningless  result,  and  it  will  be  necessary  to  constrain  the  calculation  to 
prevent  it.  With  these  restrictions,  however,  the  procedure  outlined  permits 
us  to  use  the  source  distributions  provided  by  the  method  of  Hess  and  Smith 
to  explore  the  effect  of  additions  to  given  hulls  of  any  shape. 

It  will  also  be  possible  to  calculate  the  effect  of  additions  to  the 
volume  in  a  somewhat  different  fashion.  This  is  by  distributing  source-sink 
pairs  arranged  as  dipoles  of  zero  strength  along  the  centerline  plane,  and 
then  using  the  method  of  steep  descent,  permitting  them  to  increase  but  not 
to  decrease  in  strength.  This  is  equivalent  to  permitting  Increases  in  volume, 
but  no  decreases.  Since  no  designer  has  been  known  to  decrease  the  volume 
of  his  hull  after  first  laying  it  out,  this  is  probably  an  eminently  practical 
technique . 

5.  USE  OF  A  METHOD  OF  STEEP  DESCENT  TO  FIT  A  SOURCE  DISTRIBUTION  TO  A  HULL  FORM 

As  pointed  out  earlier,  it  is  difficult  to  find  a  source  distribution  on 
a  centerline  plane  which  describes  accurately  a  given  hull  form.  However,  it 
is  certainly  possible  to  do  so  by  the  method  of  steep  descent,  although  it 
may  not  be  an  economical  method  of  calculation.  Suppose  we  have  a  hull  that 
is  described  by  a  set  of  n  "offsets",  that  is,  by  n  values  of  the  coordinate 
y  where  the  hull  is  described  by  the  equation  y  -  t  y(x,z).  Suppose  further 
that  we  have  a  total  of  N  sources  placed  on  the  centerline  plane  which  are 
assigned  strengths  m^,  and  ve  use  these  to  provide  an  approximation  to  the 
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•hope  of  Che  hull.  If  we  can  wrice  the  y-coordinate  of  the  closed  stream¬ 
line  corresponding  to  this  approximation  as  Y^,  and  the  corresponding  value 
for  the  actual  hull  as  y^,  we  way  define  a  function 


S  -  21  (Y,-y  )' 

s-i  * 


...(4.22) 


The  function  S  is  a  measure  of  the  deviation  of  the  approximation  from  the 
actual  hull.  If  we  can  find  the  relationship  hY  ,  k  •  1,...,N  for  each  of 

K 

the  Yf,  we  can  write 


h«k 


-  Zl2(Y  -y  )  hY  . 

*“k 


k*l  |M  t  }M  a 


...(4.23) 


We  may  now  follow  the  method  of  steep  descent  to  find  the  set  of  values  of 

the  which  will  minimize  S  and  so  make  the  approximation  best  fit  the  actual 

N 

hull  in  the  sense  of  least  squares.  It  is  necessary ,  of  course,  that  £  a  «  0 

k-1  k 

in  order  to  insure  closure  of  the  streamlines.  Therefore  the  descent  must  be 

N 

in  a  direction  orthogonal  to  the  gradient  of  £!  hl  .  When  we  have  minimized 

k»l 

S  we  have  found  the  best  fit  for  the  set  of  sources  chosen.  It  is  always 
possible  to  add  more  sources  if  the  fit  is  not  good  enough,  or  to  provide 
weights  to  the  terms  in  (4.22)  and  (4.23)  if  it  is  felt  that  certain  places 
on  the  hull  are  more  important  than  others  in  getting  a  good  fit. 

The  method  as  outlined  sounds  simple.  Unfortunately  the  relationship 
between  Y  and  m,  implied  by  hY  is  not  an  easy  one  to  write  down,  since  the 

.  K  x  .  — 

dmk 

streamlines  are  frequently  found  by  integrating  along  a  streamline.  However, 
the  method  is  possible  in  principle  and  may  be  worth  using  if  nothing  else 
can  be  found.  A  rough  approximation  to  ^Y  may  be  found  by  using  a  variation 

K 

of  Mlchell'c  formula  --  l.e.  assume  chat  Che  addition  of  width  to  the  hull  la 
correctly  expressed  by  equation  (3.2)  If  the  addition  Is  small  enough. 
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V*  Summary  and  Conclusions 


Although  some  knowledge  of  the  theory  of  wave-making  resistance  of  ships 
has  been  available  for  about  eighty  years  and  thousands  of  tests  have  been 
made  to  measure  this  resistance,  there  has  been  little  application  of  this 
theoretical  knowledge  to  decreasing  the  resistance  of  ships.  The  main  advances 
in  application  of  the  theory  have  come  in  the  past  ten  years,  and  all  of  these 
have  been  limited  in  their  usefulness.  The  principal  contributions  to  this 
advance  have  been  these: 

1.  The  observation  by  Inui  that  the  agreement  between  theoretical  and 
measured  resistance  could  be  greatly  Improved  if  the  hull  form  were  defined 
by  the  closed  streamlines  generated  by  its  definition  in  terms  of  sources 
and  sinks,  rather  than  by  the  simpler  approximation  used  earlier. 

2.  A  discovery  by  Inui  that  the  waves  of  some  ship  forms  could  be 
largely  canceled  by  installation  of  a  large  sphere  beneath  the  bow  and  another 
beneath  the  stern.  Wigley  investigated  this  problem  and  narrowly  missed  the 
discovery  some  twenty-five  years  earlier  [20^. 

3.  A  series  of  applications  of  the  Rite  method  by  Weinblum  and  others, 
and  an  employment  of  the  calculus  of  variations  by  Karp,  Kotik,  and  Lurye 

to  find  hull  forms  of  decreased  wave-making  resistance.  These  applications 
have  served  principally  to  dispel  the  belief  that  nothing  could  be  done. 

In  addition,  there  have  been  several  advances  in  related  fields  which 
have  not  yet  been  fully  exploited,  but  which  appear  to  have  the  potential 
to  make  a  decisive  contribution  to  the  ability  to  calculate  hull  forms  of 
decreased  wave-making  resistance.  Perhaps  the  most  Important  is  the  general 
availability  of  high-speed  computers,  which  make  possible  the  employment  of 
mathematical  methods  which  would  be  Impractical  without  them.  There  has 
also  been  developed  the  mathematical  method  of  steep  descent,  which  appears 
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to  be  Ideally  suited  for  finding  hull  forms  of  decreased  wave-making 
resistance.  An  example  is  calculated  in  part  IV  and  compared  with  earlier 

> 

results  by  the  calculus  of  variations.  The  method  of  steep  descent  can  be 
used  in  a  hand  calculation  in  simple  cases ,  as  in  the  example  cited  in  part 
IV,  but  a  high-speed  computer  is  needed  in  more  complicated  problems  of  ships 
of  a  practical  shape.  In  addition  the  ability  to  describe  an  arbitrary  hull 
form  in  terms  of  sources  and  sinks  has  been  developed  by  Smith  and  Hess  in 
the  past  three  years ,  and  the  conversion  of  the  wave-making  resistance 
equations  into  a  form  which  will  permit  exploitation  of  this  ability  with 
the  method  of  steep  descent  is  shown  in  part  II  of  this  paper.  In  consequence, 
it  appears  that  the  mathematical  machinery  necessary  to  obtain  much  more 
Improvement  than  has  been  possible  heretofore  is  now  available.  The  only 
obstacle  to  a  completely  general  method  of  finding  improved  hull  forms  is 
the  lack  of  a  simple  way  to  find  a  source  distribution  within  the  confines 
of  the  hull,  rather  than  on  its  surface,  which  generates  an  arbitrary  hull 
form.  However,  the  method  of  Smith  and  Hess  for  finding  a  source  distribu¬ 
tion  on  the  surface  of  the  hull  equivalent  to  an  arbitrary  hull  permits  the 
easy  calculation  of  any  change  which  adds  to  a  hull  volume,  and  in  view  of 
the  normal  method  used  for  designing  ships,  this  should  nearly  always  be 
sufficient.  A  method  for  carrying  out  such  a  calculation  is  outlined  in 
part  IV. 

In  summary,  it  appears  that  although  little  practical  application  has 
yet  been  made  of  mathematical  methods  for  finding  hull  forms  of  decreased 
wave-making  resistance,  the  techniques  are  now  available  which  will  permit 
this  to  be  generally  done. 
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Appendix  1 . 


Integration  of  the  Resistance  Equations* 


In  this  section  estimates  will  be  made  of  the  value  of  the  smoothly 
varying  portion  of  the  wave-making  resistance  of  an  assemblage  of  sources 
and  sinks,  and  then  a  closed  form  solution  will  be  obtained  for  this  portion 
of  the  resistance*  After  this  a  limit  will  be  derived  for  the  error  incurred 
by  terminating  the  numerical  integration  of  the  expression  for  the  fluctuat¬ 
ing  part  of  the  resistance  at  some  zero  of  the  integrand. 

The  discussion  here  will  be  based  on  equations  (2.17),  (2.18),  (2.19), 
(2.21),  and  (2.22).  For  convenience  they  are  reproduced  here: 


R  *  l6ff< 


n  n  -1 

+  2  H  Z  mf 

s»r+l  r*l 


K  m  2e"Prr  |  <»~prrt(l+t)1/2t”1/2dt  + 

0  I  r=l  r  2  o 

^rr!rras  e'p™  J5e*p”t(l+t)l/2t"1/2cos^rs(l+t)1/2j  coS^'rstl/2(l+t)1/^ 


...(2.17) 


Here  we  have  act 

p  a  2R  f  ;  p  =  K  (f  4f );q  *=  K  (h  -h  );  and  q'  ®  K  (k  -k  ) , 

rr  o  r  rrs  ox  r  sJ,Mrs  ov  r  3  *  Hr3  ov  r  s' 


R  =  +  R*2^ 


..(2.18) 


£  k 


...(2.19) 


where 


Ko2inr2c’Prr  J  e-prrt(l+t)l/r2t-l/2dt 


...(2.19a) 


n  n-l 
2  21  Z  R 
8»r+l  r=l 


!6fljpK02  J"  c~Pr=t(l+t)1/2  t_1/2cos  ^r3(l+t)1/2J 


...(2.21) 


(l+t)J 
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The  expression  rO)  Is  the  smoothly  varying  part  of  the  resistance,  while 

R^^  is  the  fluctuating  part.  Each  term  R^  of  the  fluctuating  part  will 

become  alternately  positive  and  negative  as  KQ  Increases  monotonlcally . 

Each  term  Ryr  of  the  smoothly  varying  part  R^'  will  start  small  for  small 

increase  monotonlcally  to  a  maximum,  and  then  decrease  as  K  becomes  still 

o 

larger. 

We  can  estimate  the  behavior  of  R  as  follows: 


rr 


l  .1/2  -1/2  J 

j  (1+t)  t  dt  + 
o 


S 


2  ( +  J  jm  dt  .  i/r  +  jm. 

o  3  ~ 

It  follows  that 


*3p 


rr 


rr 


» 


rr 


-  l6*f*o  mr  | 


rr 


[ 


44 


473  e 
P 


M  ■ 


-  IhrpK  V  ‘’2K°fl 
'  °  r  2 


44 


r  UJT  +  473  e 


2Kofr  - 


...(A.l) 


For  any  fr  >  0,  this  function  is  small  for  small  K0,  Increases  smoothly  to 
a  maximum t  and  then  decreases  to  aero  as  K0  becomes  infinite.  Since  our 
derivation  of  the  velocity  potential  has  assumed  that  f  >  0  so  that  the  flow 


r 

can  be  described  by  a  source  and  its  mirror  image  in  the  surface p  the  restrict 
tion  on  fy  adds  no  restriction  not  already  implicit  in  the  formulas  for  the 
resistance.  The  relation  ■  g/c  means  that  sero  speed  corresponds  to 
infinite  Rq  and  vice  versa.  The  maximum  value  of  Rrr  as  estimated  in  equation 
(A.l)  is  reached  between  the  values 

4  <  2K  fr<  6  . 

0  This  corresponds  to 


...(A.2) 
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However,  this  approximate  result  really  applies  to  the  coefficient 
of  resistsnce  defined  by  Crr  -  Rrrc"2,  rather  than  the  resistance  itself. 
This  is  because  the  definition  of  source  strength  used  here  include#  a 
factor  c  so  that  for  .a  given  hull  source  strength  increases  with  increas¬ 
ing  speed.  See  footnote  to  page  10. 


Equations  (2,19)  and  (2.19a)  show  that  t  the  smoothly  varying  por¬ 
tion  of  the  resistance,  is  for  a  given  speed  and  strength  of  a  source  a 
function  only  of  the  depth  of  the  source.  On  the  other  hand,  equations  (2.21) 
and  (2.22)  show  that  the  fluctuating  part  of  the  wave-making  resistance  is 
a  function  not  only  of  the  sum  of  the  depths  of  each  pair  of  interacting 
sources,  but  also  of  their  fore-and-aft  and  athwartships  displacement  from 
each  other*  This  arises  through  the  quantity  qra  which  is  a  function  of 
J  ^r*^s  |  *  anc*  t*ie  quantity  q£g  which  is  a  function  of  |  kr-kg  J  ,  the  fore- 
and-aft  and  athwartships  separations  respectively.  These  results  provide  a 
simple  explanation  of  a  phenomenon  apparently  first  noticed  by  Inui.*  He 
found  that  to  bring  the  calculated  wave  resistance  of  a  ship  into  conformity 
with  the  observed  wave  resistance  he  had  to  adjust  the  position  of  the  waves 
originating  at  the  stern  of  the  ship  aft  of  their  actual  position  in  the 
fluctuating  term  of  his  formula  for  wave  resistance,  but  he  needed  to  moke 
no  such  adjustment  in  the  smoothly  varying  portion.  It  appears  that  the 
smoothly  varying  portion  is  not  a  function  of  horizontal  position,  so  it  is 
insensitive  to  corrections  made  in  that  position. 

*lnui,  Takao,  "Study  on  Wave-Making  Resistance  of  Ships",  The  Society  of 
Naval  Architects  of  Japan,  60th  Anniversary  Series.  Vol.  2.  pp.  172-355.  The 
discussion  starts  on  page  207. 
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The  Integral  of  equation  (2. 19a)  may  be  evaluated  in  closed  form  as  a 


Laplace  transform.  From  Tables  of  Integral  Transforms,  Vol  I,  Bateman 
Manuscript  Project,  p,  139,  Section  4,3,  transform  (17),  we  may  find  Its 
value.  The  transform  as  tabulated  is  this: 

If  f (t)  «  0,  0  <  t  <  2b 

-  (t+2*)^  (t-2b)“^  ,  t  ->  2b, 

and 

arg(a+b)  <tt,  Re  \)  <1  1 , 


then  the  Laplace  transform  is 

LffCO]  -  g(p)  -  W-*)]  P-1  e“*a+b*p  k2l/  J(«+b)pJ  . 

We  can  put  the  integral  of  equation  (2,19a)  in  this  form  if  we  set  l/  *  1/2, 
a  *  1/2,  b*0,  arg(a+b)  ■  0,  Then,  since  csc(v/2)  ■  1,  it  follows  that 


/•“Ptrt(l+t)1/Vl/2  dt  -  (  t/2)  prr_1  e'P«/2  k^p^/2).  ...(A.3) 

o 


If  we  put  this  result  into  equation  (2,19a)  we  get 


R 


rr 


-  4?r2  pK02Prr"1 


2  _-3prr/2 


m  e 
r 


kx(prr/2) 


,..(2.20) 


The  function  k^(Prr/2)  is  Bateman's  function.  It  is  a  composite  of  known 
hypergeometrlc  functions  and  the  gamma  function,  and  can  certainly  be  tabula** 
ted  without  difficulty. 

There  does  not  appear  to  be  a  tabulated  Laplace  transform  for  the  fluc¬ 
tuating  terms  of  the  wave -resistance  R^.  However,  we  can  integrate  the 
function  numerically,  and  it  is  possible  to  show  that  if  q^  *  0  and  we 
Integrate  up  to  any  cero  of  the  factor  cos  q^^/l+'t  except  the  first  one  the 
error  is  less  than  the  value  of  the  Integral  between  that  aero  and  the  next 
one.  To  see  this,  let  1+t  -  v  In  equation  (2.22).  Then  we  have  Cor  the  Integral 
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R^  *  const 


(A.4) 


J  e"v  Rr>  (v^-1)  cos(q  v)  dv. 
1  rB 


We  aay  write  this  as  an  alternating  aeries: 


Rrs  “  const{  I  f  l  - 

r 

+ 

v- 

j3 

_ 

+ 

U1 1 

vi 

v2 

v3 

'i 

e  *>r,(v2«l)‘l/2coi(qr|v)  dv 
...(Ae5) 


The  limits  V2*  v^,  , •  •  are  the  successive  zeros  of  cos(qrsv)a  If  we  can 

show  that  the  terms  after  some  point  become  progressively  smaller  and  tend  to 

zero,  then  we  will  have  shown  not  only  that  the  Integral  of  equation  (A.4) 

and  hence  of  (2.22)  converges  but  that  the  error  made  by  dropping  all  terms 

beyond  some  point  is  less  than  the  magnitude  of  the  first  term  dropped.  This 

requires  that  the  first  term  dropped  be  smaller  in  magnitude  than  the  previous 

one.  We  observe  that  if  we  break  up  two  successive  terms  in  (A. 5)  Into  small 

Intervals  of  equal  argument  of  the  cosine,  then  not  only  will  the  intervals 

At  of  the  later  term  be  no  greater  for  equal  argument  than  the  corresponding 

Intervals  of  the  earlier  term,  but  the  absolute  value  of  the  integrand  will 

be  smaller  in  every  case  for  a  small  interval  in  the  later  term  than  for 

the  corresponding  interval  in  the  earlier  term.  From  this  it  follows  that 

every  term  of  (A. 5)  except  perhaps  the  first  one  is  certainly  followed  by  a 

term  which  is  smaller  in  Absolute  value.  But  if  this  is  so,  then  we  may 

terminate  the  Integration  of  (2.22)  when  q'  ■  0  at  any  zero  of  the  cosine, 

r  s 

except  perhaps  the  first  one,  with  the  knowledge  that  the  error  is  less  than 

the  value  of  the  integral  between  that  zero  and  the  next  one. 

If  we  set  q  *  0,  q'  ^0,  then 
nrs  nrs 

R  *  const  Je  V  (v2-l)-1^  cos(q'  v  n/v^I  )  dv  ..»(A.6) 

rs  ^  rs 

and  the  identical  argument  holds  as  for  q'  *  0,  q  4  0.  For  the  case  where 

rs  rs 

neither  qra  nor  q£#  is  zero,  the  estimate  of  the  error  to  the  precision 
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obtained  above  la  not  easy  to  obtain,  but  a  somewhat  less  precise  upper  limit 

of  the  error  is  worked  out  in  the  text* 

It  is  possible  to  draw  from  equations  (2* 19a)  to  (2*22)  some  simple  > 

deductions  which  show  that  they  provide  known  results  in  limiting  cases* 

For  example,  if  we  let  q^  *  q£fl  -  0,  which  means  two  sources  superposed, 

we  see  that  for  m  ■  hb  we  get 
r  s 

2  2  2 

R  *  R  +  2R„_  +  R.„  *  const  x  (m  -2m  +  m  )  x  integral  »  0. 

rr  rs  ss  r  r  r 

Similarly  if  we  have  two  positive  equal  sources  superposed  we  get  a 

resistance  four  times  as  great  as  for  a  single  one  --  which  would  follow 

from  the  fact  that  superposing  two  equal  sources  is  the  same  as  doubling  the 

strength  of  a  single  source,  and  the  smootnly  varying  portion  of  the  resistance 

varies  as  the  square  of  the  source  strength* 

It  19  also  possible  to  derive  from. these  equations  the  effect  of  placing 

a  source  and  an  equal  sink  near  to  and  below  a  given  source*  With  proper 

separation  and  sufficient  depth  below  the  surface  of  both  the  given  source 

and  the  added  source-sink  pair,  we  can  get  as  good  cancellation  as  we  please 

of  the  given  source's  wave  pattern*  This  follows  because  with  large  p^, 

which  results  from  large  depth  f^  of  the  given  source  and  f^  of  the  added 

sources,  very  nearly  all  the  value  of  Integrals  in  equations  (2.19a)  and 

(2,22)  is  obtained  for  very  small  t*  We  need  only  choose  q  such  that 

rs 

cos  q  *  1  and  adjust  the  value  of  m  to  get  virtually  complete  cancellation 
rs  s 

if  fr  and  fs  are  large  enough* 

It  may  be  desirable  in  some  cases  to  use  dipole  moments  rather  than 
sources  and  sinks  to  describe  the  ship*  If  this  is  done,  then  the  resistance 
equation  and  all  that  follows  from  it  can  be  placed  in  a  form  which  is  parallel 
to  that  used  for  sources  and  sinks  in  equation  (2*17)*  If  the  dipole  moment 
is  Mr,  then  the  resistance  becomes 
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R  -  !6rpK 


,  r  °  2  - 

«r  i 

I  r-1  r  2 


j  e’Pr*t(l+t)3/2t’1/2  dt  + 


+  2  £  r  e’P*8  J  e'Prst(l+t)3/2t"l/2  cos  q  (1+t)1 

s-r+1  r-1  ^ 8  2  o  JB 


H- 


.•.(A. 7) 


The  argument  concerning  the  error  incurred  bv  termination  the  integral  ia  a 
little  more  complicated  than  for  the  case  cf  3ources  and  sinks,  but  similar 


in  nature:  it  is  clear  that  there  is  some  value  of  t  beyond  which  the  error 


incurred  by  terminating  the  integral  (for  q*  -  0)  at  any  zero  of 

rs 

cos  k.U«>1/2]  is  less  than  the  value  of  the  integral  between  that  zero 
and  the  next  one. 
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